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EDITORIAL 


Las preocupaciones por el estado del medio ambiente son ahora un tema cotidiano en las 
noticias que recibimos en la prensa escrita y por los medios electrónicos, incluyendo el 
internet Esto es en sí una buena noticia pues significa que más gente se entera y se inte- 
resa por esa situación y presumiblemente, más gente habrá en posición y con la actitud 
de tomar medidas preventivas y correctivas. Es positivo que esta información se difunda 
de manera amplia pues la labor de cuidar el medio ambiente no es labor de una sola 
persona o comunidad, sino de todos. En la medida en que todos colaboremos, mayores 
probabilidades habrán de alcanzar un estadio más adecuado que el que se presenta 
actualmente. En estos dia del verano se celebra al ‘Día Mundial del Medio Ambiente’ lo 
que representa un momento propicio para reflexionar acerca de este importante tema. 

En este entorno, la labor de los Investigadores y de los editores de revistas científicas se 
magnifica y se hace más necesaria para difundir de manera amplia y con un sano criterio 
editorial, los resultados de la investigación llevada a cabo generalmente con recursos 
públicos. Son precisamente tos investigadores quienes llevan sobre sí la enorme respon- 
sabilidad de generar conciencia entre la opinión pública, de una manera inteligente y 
ordenada, para facilitar la interpretación de los efectos que sobre el medio ambiente 
tienen las prácticas humanas. Muy relacionado con este aspecto, está la situación de la 
productividad de ia tierra y la rentabilidad de los negocios derivados de las diversas acti- 
vidades de la sociedad. La investigación aspira a mejorar las condiciones de vida de la 
humanidad (excepto la investigación militar, que se afana por buscar métodos eficaces 
para aniquilar a los enemigos); en este entorno, la investigación forestal y sobre el apro- 
vechamiento de los recursos forestales en un país como México, tiene el fuerte compro- 
miso de propiciar un uso racional de uno de los mayores recursos naturales renovables 
de la humanidad: el bosque y sus productos y servicios derivados. 

En este número nos complace presentar los resultados de varias investigaciones rele- 
vantes. El primer tema está relacionado con los servicios ambientales de una cuenca 
forestal en el occidente de México. El aprovechamiento eficiente de esos servicios parece 
ser una de tas mejores apuestas que puedan hacerse en pro del medio ambiente y de tos 
seres humanos que habitan en los diversos ecosistemas. Seguramente que los datos 
registrados en este estudios darán mejores herramientas a los tomadores de decisiones 
para tratar de aprovechar estas posibilidades de manera sabia. 

El siguiente artículo es un esfuerzo Interesante por explicar un problema que enfrenta la 
producción nacional de madera aserrada con la pérdida de competitividad frente a la 
madera importada. Con el objetivo de determinar los factores que podrían revertir esta 
tendencia, se usó un modelo de programación no lineal. Los resultados indican que una 
disminución de los costos de producción, incluyendo los de transporte y la aplicación de 
un arancel reduciría las importaciones, respecto al nivel observado en el periodo 2001- 
2003, y por otra parte incrementaría la producción. Esta clase de estudios aportan datos 
útiles para los encargados de dirigir la política nacional y de hecho, para todos los intere- 
sados en el sano desarrollo del sector forestal mexicano. 

Una condición ineludible para el desarrollo forestal equilibrado y sustentabie de un país, 
es el tema de la re forestación. No se puede aspirar a alcanzar un crecimiento sano de la 
industria forestal si no se cuenta con métodos funcionales y eficientes para el crecimiento 
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de las plántulas y los nuevos árboles que reemplacen a los que se corten e incrementen, 
con el tiempo, la masa forestal de un bosque natural o plantación. Un artículo en este 
número pretende justamente eso, al caracterizar diversos medios de crecimiento para la 
producción forestal en vivero. Es un esfuerzo importante y valioso que aporta datos de 
suma utilidad para los profesionales de los viveros y las dependencias a cargo de la refo- 
restación. 

Los dos últimos manuscritos se hicieron en uno de ios estados más pobres de ia Repú- 
biica Mexicana, ubicado en el sur del país, el estado de Oaxaca. Se trata de temas de 
investigación sobre tres de las principales especies maderables del país: el Pínus patula, 
la Cedrela odorata y otra especie tropical representativa de los bosques tropicales de esta 
nación, Tabebuia donnell-smithii. Son tres especies muy apreciadas en el comercio y que 
tienen una iarga historia en México y en otras partes del mundo. Que las tres existan en 
un solo estado de la república es una indicación de la enorme biodiversidad de que disfru- 
tamos en muchas regiones de México, pero que como en muchos otros países en vías 
de desarrollo, no hemos sabido aprovechar a su máxima capacidad en un contexto 
sustentable: por el contrario, podemos hablar más bien de una situación de explotación y 
depredación. Los trabajos aquí presentados nos ilustran sobre la estructura del bosque 
bajo manejo en esas regiones, así como sobre el crecimiento y rendimiento maderable 
de las especies mencionadas. La información que contienen es de gran valor y utilidad 
para quienes tengan interés en manejar y aprovechar racionalmente sus bosques. Es con 
profunda tristeza que debemos mencionar que uno de esos trabajos, es un producto 
póstumo de uno de los investigadores forestales más destacados de esta generación en 
México y que lamentablemente falleció en octubre de 2007, el Dr. Miguel Ángel Musáiem 
Santiago. Que los resultados de su trabajo, así como de los múltiples trabajos que ante- 
cedieron a éste, sirvan a la nación, será el mejor homenaje que podamos hacer a su 
memoria. En el Comité Editorial de esta revista con ia que colaboró por muchos años, 
lo vamos a extrañar. 

Concluyo este editorial haciendo votos porque se redoblen los esfuerzos por proteger el 
medio ambiente, al mismo tiempo que se busca el máximo beneficio para ios que viven 
de este importante recurso, así como para los que dependen de sus industrias derivadas. 
Si no actuamos de manera inteligente y coordinada, el juicio de la historia no nos será 
benigno. Seguramente, no nos absolverá. 


Raymundo Dávalos Sotelo 

Editor 
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lARTICULO DE INVESTIGACION 


Oferta hídrica de la cuenca forestal 
Tapalpa, Jalisco, orientada hacia los 
servicios ambientales 

Water supply of the Tapalpa forestry watershed in Jalisco, 
focused to environmental Services provisión 

Juan de D. Benavides-Solorioi, 

Manuel de J. González-Guillén2, Cristopher López-Paniagua'' 

y J. René Valdez-Lazalde2 

RESUMEN 


Con la finalidad de aportar elementos para determinar la factibilidad de crear un mercado de 
servicios ambientales hidrológicos, se estimó el balance hidrológico en la cuenca forestal de Tapalpa, 
Jalisco, México, aplicando el modelo hidrológico SWAT. Éste presentó un buen ajuste (R2 = 0,85) con 
la producción de agua, al ser calibrado en una de las cinco subcuencas que integran el área de 
estudio. Las subcuencas forestales presentan un buen estado hidrológico al permitir la infiltración del 
agua al suelo a través del escurrimiento subsuperficial para luego emerger como escurrimiento super- 
ficial, con poca producción de sedimentos. Las subcuencas con mayor superficie agricola o pecuaria 
presentan mayor escurrimiento superficial y mayor producción de sedimentos; la evapotranspiración 
en las subcuencas forestales es también más elevada. La oferta total del recurso hídrico, para un año 
promedio, asciende a 42,963,900 m3 año-1. 

PALABRAS CLAVE: 

Balance hidrológico, cuenca forestal, modelo swAT, servicios ambientales hidrológicos. 


ABSTRACT 


In order to assess the feasibility of creating a water market in the Tapalpa, Jalisco, México 
watershed, the SWAT model was fitted to obtain the watershed water balance. The model was cali- 
brated in one of the five subwatersheds included in the Tapalpa watershed with a good fit (R2 = 0,85). 
All five forestry subwatersheds had good hydrologic behavior, allowing high infiltration water rates, 
which later on emerges as surface runoff with low sediment yieids. The subwatersheds with high 
percentage of agricultura and grassiands yieid more surface runoff as well as more sediment yieids. 
The evapotranspiration from the forestry subwatersheds were higher than agricultural subwatersheds. 
The water supply from the Tapalpa watershed is 42’963,900 m3 per year on average. 

KEY WORDS: 

Water balance, forestry watershed, SWAT model, hydrologic environmental Services. 


1 INIFAP-CIRPAC, Parque Los Colomos s/n. Col. Providencia, C.P. 44660, Guadalajara, Jalisco, México. Tel. 
(33)4641-2061 Ext. 118, correo: benavides.juandedios@inifap.gob.mx, cristopherlp@colpos.mx. 

2 Postgrado Forestal. Colegio de Postgraduados. Km. 36.5 carr. México- Texcoco. Montecillo, Estado de 
México. 56230. correo: manuelg@colpos.mx, valdez@colpos.mx. 
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INTRODUCCIÓN 

La distribución irregular de la 
cantidad y calidad de agua dulce, aunada 
a problemas de contaminación de la 
misma, representan un problema que se 
ha acentuado en los últimos años en 
varias regiones de México. En este 
sentido, existen evidencias alarmantes de 
sobreexplotación de acuíferos y contami- 
nación de cuerpos de agua y escurri- 
mientos superficiales. Informes de la 
Comisión Nacional del Agua (cna) indican 
que dos tercios de los 188 acuíferos de 
mayor importancia del país están siendo 
sobreexplotados (dof, 2003), debido a 
que el promedio de extracción anual de 
agua supera aproximadamente en dos 
veces el nivel de recarga natural. 

Un ejemplo de problemas de calidad 
de agua superficial lo constituye la cuenca 
hidrológica Tapalpa, localizada en el estado 
de Jalisco. Este deterioro puede ser debido 
al aumento poblacional en la región y a un 
cambio de uso del suelo distinto al forestal, 
propiciando más áreas desprovistas de 
vegetación. Además, las características 
propias de la región han permitido un 
desarrollo en la actividad turística a tal 
grado que empiezan a vislumbrarse 
problemas en la población de Tapalpa, prin- 
cipalmente relacionados con la provisión de 
servicios básicos -entre ellos el recurso 
hídrico- en el cual se exhibe una falta de 
planeación del manejo del agua. 

Lo anterior deja de manifiesto la 
necesidad de implementar estrategias 
tendientes a mejorar el aprovechamiento 
y conservación del recurso hídrico 
mediante la conservación de la cubierta 
forestal en su estado natural; cumpliendo 
la función de regulador hídrico (Chang, 
2003) al permitir altas tasas de infiltra- 
ción, y en consecuencia a la recarga de 
mantos freáticos con agua de alta 
calidad. Una de estas estrategias en 
México, desarrollada por el Gobierno 


Oferta h id rica de la cuenca forestal Tapalpa... 

Federal a través de la conafor, la consti- 
tuye la creación de mercados de servicios 
ambientales hídricos (sah), a través de 
los cuales se busca la retribución o 
compensación monetaria por parte de los 
usuarios de los sah a los dueños de los 
ecosistemas que los proporcionan. A nivel 
mundial, los mercados de sah son un 
mecanismo creciente para conservar y 
rehabilitar los recursos naturales provee- 
dores de beneficios públicos (Torres et 
al., 2005). Un programa pionero de este 
tipo se desarrolló en Costa Rica en 1997 
(Pagiola, 2003); otros programas pilotos 
se han desarrollado en Latinoamérica y 
Asia (Echavarria, 2003; Kerr, 2003). En 
años recientes, estos programas se han 
introducido en países en desarrollo con 
resultados exitosos (World Bank, 2005). 

En México, desde hace varias 
décadas se han venido desarrollando una 
serie de iniciativas sobre negociaciones e 
intercambios entre usuarios y proveedores 
de servicios ambientales con el objetivo de 
garantizar o mejorar la provisión de tales 
servicios. Estos intercambios no eran 
conceptualizados como “mecanismos de 
pago de sah” o como un mercado formal, 
sino como simples arreglos entre partes 
interesadas (Madrid, 2005). 

El éxito de un mercado de sah se 
basa en la creación de un sistema en el 
cual interactúan de manera armónica las 
fuerzas de oferta y demanda, lo cual se 
logra al contar con información que 
permita su correcta determinación y 
caracterización. Por el lado de la oferta 
hídrica, el estudio de los aspectos físicos 
en la zona de estudio, incluyendo el cono- 
cimiento de las salidas de la misma (p. ej. 
el balance hídrico y sus componentes) es 
una parte que coadyuvará al éxito en la 
implementación de un mercado de sah 
(Pagiola et al., 2003). 

Un estudio de la oferta de sah debe 
considerar el balance hídrico de las 
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cuencas, conocimiento que se obtiene 
mediante la toma de datos de campo y el 
uso de modelos hidrológicos. Estos 
últimos permiten conocer algunas varia- 
bles que no hayan sido tomadas en el 
campo y que son fundamentales para 
realizar un balance hidrológico completo y 
confiable. Los modelos hidrológicos 
tienen varias ventajas con respecto a los 
métodos de medición directos o de 
campo, tales como: realizar proyecciones 
futuras con base en los cambios de uso 
del suelo minimizando el impacto sobre 
los ecosistemas, ajustar variables de 
medición de campo, además de la 
ventaja de la economía, ya que permiten 
tomar datos para calibrar modelos y 
obtener resultados sin esperar largos 
periodos de tiempo en la toma de datos 
de campo (Brooks et al., 2003). 

Los modelos empleados para la 
simulación de cuencas son representa- 
ciones simplificadas de los sistemas hidro- 
lógicos reales que permiten estudiar el 
funcionamiento de los mismos y sus 
respuestas a los diferentes factores que se 
incorporan al modelo. Además predicen la 
respuesta hidrológica ante modificaciones 
de sus componentes, cuando el modelo 
está bien sustentado. La mayoría de estos 
modelos se basan en aproximaciones de 
sistemas y difieren en la forma de cómo se 
considera cada componente del ciclo 
hidrológico (Brooks et al., 2003). En 
general, los modelos hidrológicos de 
cuencas tienen cinco componentes 
básicos: (1) los procesos y características 
hidrológicas de las cuencas; (2) los datos 
de entrada: (3) las ecuaciones gober- 
nantes: (4) los límites y condiciones 
iniciales, y (5) las salidas (Singh, 1995). 

El modelo swat, utilizado en este 
estudio, es un modelo matemático, física- 
mente basado y con propiedades deter- 
minísticas; se considera un modelo 
continuo, no diseñado para predecir 
efectos por tormentas individuales. Tuvo 


sus inicios en el Servicio de Investigación 
Agrícola (ars) que forma parte del Depar- 
tamento de Agricultura de los Estados 
Unidos (usda) cuyo principal promotor fue 
Jeff Arnoid en la década de 1990. El 
modelo antecesor al swat fue el modelo 
hidrológico swrrb (Simulator for Water 
Resources in Rural Basíns) desarrollado 
por Jim Williams y colaboradores 
(Williams etal., 1985; Arnoid etal., 1990). 
El SWAT fue diseñado para predecir el 
impacto de algunas prácticas de manejo 
aplicadas al suelo en la producción de 
agua y de sedimentos en cuencas 
grandes y complejas. Integra variables 
del clima, suelo, topografía, vegetación y 
prácticas de manejo del suelo, para 
predecir los cambios que pueden ocurrir 
en una cuenca. 

Debido a su versatilidad y eficiencia, 
el modelo swat ha sido ampliamente utili- 
zado en varios lugares del mundo como 
apoyo en la toma de decisiones sobre el 
manejo del agua (Maldonado de León et 
al., 2001; Hernández, 2003; Benavides- 
Solorio et al., 2005a). Se han realizado 
balances hidrológicos en Italia, Alemania, 
Australia, Estados Unidos, Sudamérica y 
México. En la mayoría de los casos el 
modelo se ha ajustado bien a los datos, 
por lo que es una buena alternativa para 
representar ciclos hidrológicos y predecir 
la cantidad de agua superficial y el 
arrastre de sedimentos que se puede 
esperar en una cuenca (Benavides- 
Solorio et al., 2003). 

En México, versiones anteriores al 
modelo swat (p. ej. el modelo swrrb), se 
han aplicado por Fernández et al. (1996) y 
Benavides-Solorio (1998) en la cuenca 
agrícola “El Tejocote” en el Estado de 
México y en la cuenca forestal "El Carrizal” 
en Jalisco, respectivamente. El primer 
autor obtuvo coeficientes de determinación 
de 0,80; 0,88 y 0,61 para tres años conse- 
cutivos en producción de agua y coefi- 
cientes de 0,79; 0,74 y 0,41 para la 
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producción de sedimentos en los mismos 
tres años consecutivos (Fernández, 1996; 
Fernández et al., 1999). Benavides-Solorio 
(1998) obtuvo coeficientes de determina- 
ción de 0,91 para la producción de agua y 
coeficientes de 0,78 para la producción de 
sedimentos en dos años consecutivos. 
Recientemente, con el uso del modelo 
SWAT se han obtenido buenas experiencias 
como la mostrada por Torres-Benites et al. 
(2004), en la cuenca “El Tejocote” con 
coeficientes de determinación de 0,94 y 
0,84 para la producción de agua. 

OBJETIVO 

Con la finalidad de aportar elementos 
para determinar la factibilidad de crear un 
mercado de sah, el presente estudio tuvo 
como objetivo cuantificar y caracterizar el 
balance hídrico total y el agua disponible 
para uso (oferta), en la cuenca Tapalpa, 
Jalisco. En un articulo complementario a 
este trabajo se reporta información refe- 
rente a la demanda hidrológica en la cuenca 
(López-Paniagua et al., 2006); información 
igualmente importante para la definición de 
un estudio de mercado hidrológico. 

METODOLOGÍA 

Descripción del área de estudio 

La cuenca Tapalpa está ubicada en 
el estado de Jalisco en la región centro- 
suroeste en el municipio del mismo 
nombre, a 130 km de la ciudad de Guada- 
lajara (Figura 1). Sus coordenadas 
extremas son: 2'1 99,500 y 2'21 7,500 L. 
N. y 621,500 y 639,500 L. O. La cuenca 
se encuentra en el sistema montañoso 
conocido como “Sierra de Tapalpa”, 
dentro del Eje Neovolcánico. Cubre un 
área de alrededor de 210 km2; 44% de su 
superficie es de uso forestal, seguido por 
28% de uso agrícola y urbano, y 27% de 
uso pecuario. 
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Perfil sociodemográfico 

De un total de 15 480 habitantes del 
municipio de Tapalpa, 7 901 (51%) son 
mujeres y 7 579 (49%) son hombres 
(iNEGi, 2001). La tasa media anual de 
crecimiento de la población pasó de 
1,76% en la década de los 70 a 2,5% en 
los 90. El número de habitantes creció en 
un 9,8% de 1995 al 2000 (emm, 2000). 

El porcentaje de la población 
urbana y rural del municipio no se ha 
modificado sustancialmente; en la 
década de 1990, 63% de la población 
era rural y el 37% urbana, para el censo 
de 1995 la población rural descendió a 
60% y la urbana aumentó a 40%. Para el 
año 2000 la población rural volvió a los 
valores de 1990 con 64% de población 
rural y 36% de población urbana. Esto 
último puede tener su origen en la inmi- 
gración de turistas que tienen sus casas 
de campo en los bosques de la región. 

Referente al alcance que tiene la 
población a los servicios básicos, 80% de 
los habitantes disponen de agua potable; 
68% de alcantarillado y 86% de energía 
eléctrica (emm, 2000). 

Planificación del proceso y recorridos 
de campo 

Se realizaron recorridos de campo 
en la zona de estudio con el fin de tener 
un panorama general de la problemática 
de SAH, así como de las características 
socioeconómicas de los diversos 
sectores de la región. En las visitas se 
contactó a personal de las instituciones 
ligadas al manejo del recurso hídrico, así 
como a las autoridades municipales, la 
asociación ganadera local, la unidad de 
riego “El Nogal”, y a los expertos que 
conformarían el equipo de trabajo. 
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Cuenca Tapalpa 


Figura 1. Localización de la cuenca Tapalpa en el estado de Jalisco. 


Determinación del balance hidrológico 
(oferta de agua) 

La cuenca Tapalpa se tomó como 
base para estimar la oferta hídrica, ya que 
cuenta con varios aspectos interesantes 
como son: contiene una subcuenca 
forestal donde se recolectaron datos de 
precipitación y escurrimiento que fueron 
utilizados para calibrar el modelo hidroló- 
gico; su sistema hidrológico incluye a la 
población de Tapalpa, la cual se abastece 
con el agua producida en la cuenca 
Tapalpa y se tienen varios manantiales 
que permiten entender el movimiento del 
agua en la cuenca. 

Modelos de elevación digital y fisio- 
grafía de la cuenca Tapalpa 

El modelo de elevación digital se 
generó a partir de datos vectoriales dispo- 
nibles en el Instituto Nacional de Estadís- 
tica, Geografía e Informática (inegi) a 
partir de las cartas digitales números 
E13B14 y F13D84 (INEGI, 2000). 


La información de campo y de cartas 
topográficas permite apreciar que en la 
sierra de Tapalpa se distinguen varias 
topoformas. En las partes altas se loca- 
lizan elevaciones escarpadas que 
alcanzan altitudes de hasta 2 660 m 
(Figura 2); entre ellas destacan los cerros 
El Estepeme, El Mescalillo, El Divisadero 
y El Fraile. En los lomeríos, la pendiente 
disminuye y sobre ellos se realizan activi- 
dades agropecuarias, destacándose el 
cultivo de maíz, los cereales y el cultivo 
de la papa. 

Clima. El clima de Tapalpa corresponde a 
un clima semi-templado, que tiene una 
precipitación invernal promedio de 84 mm, 
y representa entre 5% y 10% de la precipi- 
tación total. La precipitación anual 
promedio es de 890 mm, tomada del año 
1961 al año 2000 (Figura 3). La época de 
lluvias generalmente se presenta de junio 
a octubre y el número de días al año con 
lluvia es de 96. La evaporación anual es de 
1 577 mm (Ruiz etal., 2003). Los extremos 
de temperaturas están en 36 °C y -6 °C. 
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Leyenda 

Red hidrológica 


Valores de elevación 

I I 2,190 - 2,290 

i 2,290 - 2,390 
[ 2,390 - 2,400 
I 2,490 - 2,660 


Fuente: Elaboración propia con datos del inegi (inegi, 2000) y campo. 

Figura 2. Intervalo altitudinal y red hidrológica de la cuenca Tapalpa. 


Precipitación — ♦ — Temperatura 



o 

o 

re 

3 

+-* 

re 

<u 

a 

E 

0) 


0 


Fuente: Elaboración propia con datos de la estación meteorológica Tapalpa de la cna. 


Figura 3. Precipitación y temperatura media mensual (1961-2000) 

de la Estación Tapalpa. 
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Los datos del clima se obtuvieron principal- 
mente de dos fuentes: (i) la estación mete- 
orológica administrada por la cna; y (ii) de 
datos de estaciones meteorológicas portá- 
tiles instaladas en la cuenca El Carrizal 
(Benavides-Solorio, 1998). 


Toma de datos en la cuenca Tapalpa 

Estaciones meteorológicas portátiles de 
El Carrizal. Se instalaron en la subcuenca 
El Carrizal cuatro pluviómetros automá- 
ticos a inicios del año 1994 y estuvieron 
en operación durante 1994 y 1995. Se 
instaló un pluviómetro en la parte alta, 
dos en la parte media y uno en la baja de 
la subcuenca El Carrizal (Benavides- 
Solorio, 1998). Los pluviómetros son del 
tipo de cono basculante, almacenan la 
información a intervalos de 1 hasta 30 
minutos, son de la marca sigma y tienen 
una precisión de hasta 0,01" (0,25 mm). 
Fueron construidos con los estándares 
del Servicio Nacional del Clima de los 
Estados Unidos. Se programaron para 
registrar datos cada 10 minutos, por lo 
que fue posible registrar la intensidad de 
la lluvia de todas las tormentas. La infor- 
mación de lluvia se procesó por tormenta 
y por lluvia diaria. En la parte baja se 
instaló una minicomputadora para alma- 
cenar información sobre temperatura, 
humedad relativa, precipitación y altura 
del tirante de agua. Los datos registrados 
durante los dos años fueron utilizados 
para realizar la calibración del modelo 
SWAT a la cuenca El Carrizal. 

Generación de estadisticos climáticos. 
Se utilizaron los modelos wxparm y 
wxGEN de EPic (Williams, 1995) para 
obtener probabilidades de lluvia y esta- 
dísticos de temperatura y precipitación. 
Se utilizaron los datos de 1961 al 2003 de 
la estación meteorológica Tapalpa para 
generar los promedios y probabilidades 
estadísticas. Se utilizó también el modelo 
CLiMGEN para generar valores mensuales 


de radiación solar mensual, humedad 
relativa y velocidad del viento, utilizando 
los datos anteriormente mencionados. El 
generador climático del modelo swat se 
utilizó con fines comparativos y única- 
mente para estimar la probabilidad de 
lluvia. 

Hidrologia. La Sierra de Tapalpa se ubica 
en la región hidrológica 16 (rh- 16), que 
corresponde a los ríos Armería-Coahua- 
yana (inegi, 2001b). El río Armería lo 
constituyen los ríos Ayuquíla y Tuxca- 
cuesco. El principal afluente del río 
Tuxcacuesco es el río Jiquilpan, cuyo 
principal aportador es el río Tapalpa. Este 
río es el abastecedor de la presa El 
Nogal. La cuenca que abastece la presa 
El Nogal, además del río Tapalpa tiene a 
los arroyos El Carrizal, Las Ánimas, El 
Zopilote y otros más pequeños. Dentro de 
éstos, existe una gran cantidad de 
manantiales los cuales determinan que la 
mayoría de los arroyos tengan corrientes 
perennes; cuando la sequía es extrema y 
el año anterior seco se presentan como 
corrientes efímeras. Los arroyos aquí 
mencionados fueron utilizados como 
base para realizar la división de la cuenca 
por subcuencas principales. 

Medición de escurrimiento en la 
subcuenca El Carrizal. Se construyó una 
obra de aforo en la parte media de la 
subcuenca El Carrizal, aproximadamente 
en el límite entre el área forestal y el área 
agropecuaria, con la finalidad de medir 
los escurrimientos de la parte alta y 
forestal de la cuenca. La obra es del tipo 
precalibrada y de forma triangular. El 
ángulo de descarga es de 120° y tiene 
una altura máxima del tirante de 1 m, 
mismo que permite medir descargas 
máximas de hasta 2,3 m3-s-i. La caracte- 
rística del vertedor triangular es que 
permite medir con alta precisión 
descargas pequeñas y grandes. Se utilizó 
la siguiente ecuación para estimar el 
volumen de descarga (Hertzier, 1938): 
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Donde: Q = descarga medida en pies 
cúbicos por segundo: y H = altura del 
tirante del agua en pies. 

La altura del tirante de agua se 
midió de manera automática con un 
sensor de presión de fluidos (pressure 
transducer). El sensor se calibró en labo- 
ratorio hasta obtener valores de altura 
del agua comparables a los observados 
con una alta correlación (r = 0,98). Se 
utilizó la computadora Micrologger 21X, 
Campbell Scientific Inc., para registrar y 
almacenar los niveles de agua a inter- 
valos de 15 minutos. Esta computadora 
también registra información de precipi- 
tación y temperatura. La información se 
almacenó en módulos de memoria y se 
retiró del equipo en intervalos de 15 días 
a un mes. La información de escurri- 
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mientes corresponde a los años 1994 y 
1995 y se utilizó para efectos de calibra- 
ción. 

Uso del suelo 

Se utilizó una imagen de satélite 
Landsat 7 etm+ (Enhanced Thematic 
Mapper Plus) del año 2002 para generar 
información actualizada de los usos del 
suelo de la cuenca Tapalpa. Este trabajo 
incluyó etapas de campo y de gabinete. 
En una primera etapa se recorrió la 
cuenca para identificar los usos del suelo 
y establecer geográficamente puntos de 
control. Después se utilizó la imagen 
corregida geométricamente en er-mapper 
6,1 utilizando el modelo de elevación de 
INEGI (iNEGi, 2000). Posteriormente se 
realizó una clasificación no supervisada, 
utilizando la técnica denominada isodata 
(Iterativa Self-Organized Analysis Tech- 
nique) (Jensen, 1996). 



A 


Leyenda 

Vegetación 

^ Sin Vegetación 
^ Cultivos 
^ Pastizal 
^ Pastizal y arbustos 
^ Pastizal y arbolado 
I Hojosas 
^ Pino baja densidad 
I Pino alta densidad 
^ Bosque M ixto 
^ Cuerpo de agua 
1:100,000 


Fuente: Flores et al. (2005). 


Figura 4. Uso del suelo en la cuenca Tapalpa. 
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La información resultante se agrupó 
en 10 categorías (Figura 4) (i) Bosque 
Mixto (Pino y latifoliadas -hojosas-); (ii) 
Pino baja densidad; (iii) Pino alta 
densidad; (iv) Latifoliadas; (v) Pastizal 
con arbolado; (vi) Pastizal con arbustos; 
(vii) Pastizal; (viii) Áreas de cultivo; 

(ix) Área urbana o sin vegetación y 

(x) Cuerpos de agua (Flores et al., 2005). 

De las 21 000 ha que tiene la cuenca, 
el principal uso del suelo es forestal ya que 
un 44% está cubierto de bosque; en su 
mayor parte con bosque de pino puro {Pinus 
michoacana, R oocarpa, R lumholtzii, R 
leiophylla y R douglasiana), con algunas 
mezclas de encino {Quercus crassipes, Q. 
obtusata, Q. candicans, Q. castanea y Q. 
scytophylla) u otras latifoliadas {Arbutus 
xalapensis, Al ñus sp., Persea sp. y 
Crataegus sp.) en un pequeño porcentaje 
(5%). Los otros usos del suelo son el 


pastizal y el agrícola-urbano con 27% y 28% 
respectivamente, y los cuerpos de agua 
ocupan un poco más de 1% de la superficie 
de la cuenca. 

Suelos de la cuenca Tapalpa. De 
acuerdo con la clasificación de fao (1988), 
los suelos se han ordenado en cinco 
unidades (Gómez-Tagle y Chávez, 1987; 
CETENAL, 1976): andosol húmico, cambisol 
crómico, luvisol crómico, feozem háplico y 
regosol éutrico (Figura 5). 

Bases teóricas del modelo SWAT 

Balance hidrológico. El modelo swat 
permite simular diferentes procesos 
físicos de una cuenca, los cuales pueden 
ser divididos por subcuencas. El principal 
componente que simula el modelo es el 
balance hidrológico, a través del cual 
posteriormente se puede predecir el 


6225001 
l_ 


6300001 
l_ 


6375001 
l_ 


6450001 
l_ 



A 


Leyenda 

Suelos 

Cambisol 
I Feozem 
Luvisol 
I Regosol 
I Andosol 


1 : 100,000 


Fuente: Elaboración propia con datos de cetenal 1976 y Gómez-Tagle y Chávez-Huerta (1987). 


Figura 5. Unidades de suelo de la cuenca Tapalpa. 


14 

movimiento de nutrientes, sedimentos y 
pesticidas. 

El balance hidrológico se conoce 
mediante la ecuación del balance de 
agua, definida de la siguiente manera: 

SW¡ = SWq + ) 

( 2 ) 

Donde: SWj = contenido final de agua 
(mm): SWq = contenido inicial de agua en el 
día i (mm); t= tiempo (días); = precipi- 
tación en el día /; Qsup = escurrimiento 
superficial en el día i\ Ea = evapotranspira- 
ción en el día /; Wsubt - percolación en el día 
Qsubs - escurrimiento subsuperficial que 
regresa al cauce en el día i. 

Clima. Las variables del clima dentro del 
modelo son importantes porque proveen las 
entradas de humedad y energía necesarias 
para controlar el balance hidrológico. Las 
variables más importantes son la precipita- 
ción, la temperatura y la radiación solar. 

Hidrología. El escurrimiento superficial se 
obtiene cuando la cantidad de agua que 
llega al suelo sobrepasa las tasas de infil- 
tración. En suelos secos la infiltración es 
alta, pero en suelos húmedos la infiltración 
empieza a decrecer muy rápido con el 
tiempo. Cuando el suelo está completa- 
mente saturado es cuando el escurri- 
miento superficial se presenta. El escurri- 
miento superficial se calcula por dos 
métodos: la curva numérica (scs, 1972) y 
el método de infiltración (Oreen y Ampt, 
1911). El primer método utiliza la ecuación 
de la curva numérica y tablas que propor- 
cionan un valor de escurrimiento para dife- 
rentes tipos de vegetación y de suelo. La 
humedad del suelo es importante en el 
empleo de este método. El segundo 
método se basa en una ecuación de las 
tasas de infiltración que dependen de la 
conductividad hidráulica. Cuando la inten- 
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sidad de la lluvia es menor que la tasa de 
infiltración toda el agua se infiltra. La poro- 
sidad de los suelos es importante en el 
empleo de este método. 

La evapotranspiración (et) es un 
término que incluye varios procesos en 
los cuales el agua de lluvia se convierte 
en vapor. Este término incluye la evapora- 
ción desde el follaje de las plantas, la 
transpiración, la sublimación y la evapora- 
ción desde el suelo. La et es el primer 
mecanismo por el cual el agua se regresa 
a la atmósfera. Estimaciones bastante 
gruesas indican que aproximadamente el 
62% de la precipitación que llega a la 
tierra es evapotranspirada. 

Existen tres métodos que utilizan el 
modelo para conocer la et potencial: 
(1) el método Penman-Monteith que requiere 
datos de radiación solar, temperatura, 
humedad relativa y velocidad del viento; (2) 
el método Priestiey-Taylor que requiere de 
datos de radiación solar, temperatura y 
humedad relativa; y (3) el método Hargre- 
aves, el cual requiere sólo de temperatura. 

La percolación se calcula para cada 
capa de suelo; el agua se mueve cuando 
la cantidad de agua sobrepasa la capa- 
cidad de campo de cada capa. Se utilizan 
varias ecuaciones que mueven el agua 
de acuerdo a la textura del suelo y a la 
conductividad hidráulica. 

El flujo subsuperficial es aquel que 
regresa a la superficie pero primero se 
introduce en el suelo. El modelo también 
utiliza varias ecuaciones para obtener 
este tipo de escurrimiento. 

Erosión. El sedimento que predice el modelo 
se basa en la ecuación universal de suelo 
modificada (Williams, 1995), dicha ecuación 
ya enunciada fue ajustada para así poder 
utilizarla en cuencas; se basa en el escurri- 
miento superficial y el escurrimiento 
máximo, además de otras variables como la 
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cobertura, suelo y pendiente. La ecuación 
utilizada fue la siguiente: 

sed = 11,8 (Qsup * Qpico * areafjruP-56 * K^sie 
* Cusle * Pusle * (3) 

Donde: sed = es la producción de sedi- 
mento en un día (toneladas): Q^^p = escurri- 
miento superficial (mm); = escurri- 

miento pico (m3 s-i); areaf,m = área de las 
unidades homogéneas de respuesta hidro- 
lógica (ha); Kusie = factor de erosionabilidad 
del suelo; = factor de cobertura de la 
ecuación universal de pérdida de suelo; 
Pusle - factor de prácticas de manejo de la 
ecuación universal de pérdida de suelo; y 
PSusie - factor topográfico de la ecuación 
universal de pérdida de suelo. 

El sedimento también se calcula por el 
flujo lateral y por el escurrimiento subsuper- 
ficial. El modelo swat, entre otras cosas, 
también predice los rendimientos de los 
cultivos y el movimiento de los nutrimentos, 
pesticidas y la calidad del agua; sin 
embargo, estas variables no fueron conside- 
radas en el modelo para este trabajo. 

Procedimiento del uso del modelo 
swAT (interfase Arc-View) 

Inicialmente se crearon todas las 
bases de datos que utiliza el modelo, 
siendo éstas las siguientes: 

1} Mapas: Modelo de elevación digital, 
mapa de suelos, mapa de uso del suelo, 
máscara con la delimitación de la cuenca, 
mapa de hidrología con la delimitación de 
los arroyos; 

2) Bases de datos para la calibración: 
Características físicas y químicas de los 
suelos, precipitación diaria, temperatura 
máxima y mínima diaria; y 

3) Bases de datos para generación de 
información climática: Estadísticas climá- 


ticas de precipitación, temperatura, radia- 
ción solar y humedad relativa. 

Delimitación de cuenca y subcuencas. El 
modelo swat inicia con una delimitación 
automática de la cuenca, para lo cual se 
requieren el modelo de elevación digital, 
la máscara de la cuenca y un mapa hidro- 
lógico digitalizado. Este último es 
opcional pero se recomienda su utiliza- 
ción ya que realiza una mejor digitaliza- 
ción de los arroyos y de la cuenca 
completa (Figura 2). 

Con esta información se generaron 
los parámetros generales para la totalidad 
de la cuenca y por cada subcuenca: 
Tapalpa, El Zopilote, El Carrizal (parte alta 
utilizada para calibración). El Carrizal (parte 
baja) y las Ánimas.EI Carrizal se dividió en 
dos debido a que la instrumentación para 
su calibración se encuentra en la parte alta 
y el vertedor en la parte media (Figura 6). 

Uso del suelo. Se selecciona el comando 
del SWAT que involucra los tipos y usos de 
suelo y se solicita el mapa con el uso del 
suelo previamente realizado. De ahí, se 
obtienen en cuadros los usos del suelo 
contenidos en el mapa. El modelo swat 
los coloca en una pantalla y solicita que 
se les otorgue un nombre de cobertura 
que utiliza en su base de datos. Aquí es 
donde se reclasificaron algunos tipos de 
vegetación que fueron identificados en 
campo en el apartado de uso del suelo. 

Los usos definidos con el modelo 
fueron las coberturas de: agricultura, 
pastizales, praderas con gramíneas, 
praderas con arbustos, encino, pino y 
bosque mezclado y cuerpos de agua. 

Integración del uso del suelo con el tipo 
de suelo 

Este procedimiento se utiliza para 
realizar una reclasificación de los suelos 
con sus usos, y con ello se obtienen las 
hectáreas de cada agrupación. 
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622500 ™™ 630000 ™™ 637500 ™™ 645000 ™™ 



Figura 6. División de subcuencas generadas por el modelo swat 


Unidades de respuesta hidrológica (hru) 
Después de la integración de usos del 
suelo y tipos de suelo se realiza una distribu- 
ción de las unidades hidrológicas que 
incluyan uno o más suelos por cada unidad. 
Este paso permite unificar criterios por uso 
del suelo y cobertura. El escurrimiento y la 
evapotranspi ración se calculan por sepa- 
rado para cada unidad. Ello permite predecir 
con certeza los valores del balance hídrico. 
Aquí el usuario puede definir si una unidad 
hidrológica se asigna a una subcuenca, o 
por el contrario varias unidades se asignan a 
una subcuenca. Para la cuenca Tapalpa se 
eligió la segunda opción. 

Elección de datos climáticos 

El siguiente menú permite elegir si se 
obtienen datos climáticos de campo o si los 
genera el modelo swat. Si se incluyen 
datos medidos en campo se deben 


presentar los archivos de precipitación, 
temperatura, radiación solar, velocidad de 
viento y humedad relativa. Para la cuenca 
Tapalpa se introdujeron los datos diarios de 
precipitación y temperatura. Con los esta- 
dísticos del generador climático mensual 
se obtuvieron los datos de radiación solar y 
humedad relativa. 

Integración de bases de datos 

Este comando permite revisar y editar 
las bases de datos necesarias para correr 
el modelo swat. Ese modelo analiza la 
correcta inclusión de los datos y verifica si 
existen inconsistencias en ellos; cuando los 
datos son ilógicos o están fuera de rango, 
el procedimiento se detiene, indicando en 
dónde se encuentra el error. También 
define espacialmente en dónde se encuen- 
tran las estaciones meteorológicas con 
entradas de precipitación, temperatura y 
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los generadores climáticos. El modelo 
examina que todas las bases de datos 
estén consideradas y que se tenga un 
orden lógico de las variables ingresadas. 

Corrida del modelo 

Cuando se ejecuta el modelo, se 
deben incluir variables como el periodo 
de la simulación, indicando el inicio y el 
término del mismo. Para el caso de la 
cuenca Tapalpa y para realizar la calibra- 
ción del modelo, se eligió el 1 de enero de 
1994 al 31 de diciembre de 1995, periodo 
que corresponde a la toma de datos. El 
método utilizado para conocer la evapo- 
transpiración fue el “Priestiey-Taylor” ya 
que los otros métodos utilizados por el 
modelo swat como Penman-Monteith y 
Hargreaves presentaron menor ajuste a 
los datos. Se eligió la salida de datos 
mensual. 

Calibración 

Al correr el modelo se pueden 
observar los resultados desde una serie 
de ventanas en pantalla de reportes. Los 
datos se comparan con los medidos en el 
campo, se observa cuál es la diferencia y 


se plantean posibles soluciones. Se 
cuenta también con un comando de cali- 
bración en el cual algunos valores de alta 
sensibilidad pueden modificarse hacia 
arriba o abajo del valor inicial asignado. Al 
hacer las modificaciones necesarias se 
vuelve a correr el modelo las veces que 
sea necesario hasta llegar a un punto de 
calibración en que los resultados del 
modelo sean similares a los medidos en 
el campo y se obtengan resultados satis- 
factorios para el usuario. Se debe tener 
especial cuidado en que los parámetros 
modificados continúen cumpliendo con 
las características físicas. 


RESULTADOS Y DISCUSIÓN 
Estimación de la oferta hídrica 

Características generales de la cuenca y 
subcuencas 

Las subcuencas fueron delimitadas 
por el modelo swat considerando los 
arroyos más importantes. La información 
correspondiente a las subcuencas se 
muestra en la tabla 1 . 


Tabla 1. Caracterización de las subcuencas de la cuenca Tapalpa. 


Subcuenca 


Altitud (m) 


Superficie 

(km2) 


Uso DEL SUELO (%) 


Mínima 

Máxima 

media 

Bosque 

Ganadero Agrícola 

Cuerpos 
de agua 

Tapalpa 

2 000 

2 660 

2 224 

121,94 

44 

28 

26 

2 

El Zopilote 

2 040 

2 620 

2 288 

57,60 

36 

28 

36 

- 

Las Ánimas 

2 020 

2 567 

2 182 

16,52 

59 

22 

18 

1 

El Carrizal^ 

2 020 

2 420 

2 251 

7,69 

85,1 

14,8 

0,1 


El Carriza|2 









(parte baia) 

2 020 

2 270 

2 070 

5,14 

15 

36 

49 

- 


1 Esta subcuenca fue utilizada para calibrar el modelo, ya que en esta área se encuentran los instrumentos de 

medición. 

2 Esta subcuenca fue utilizada sólo para fines comparativos ya que la parte alta fue separada por la obra aforadora. 
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Calibración del modelo 

El resultado de la primera salida del 
modelo varió considerablemente respecto 
a la información medida en la cuenca El 
Carrizal: la producción de agua se sobre- 
estimó en los dos años en que se evaluó 
el modelo. Se utilizaron solamente dos 
años porque son los de duración del 
proyecto y en los cuales se tienen datos. 
Los picos máximos mensuales fueron 
sobreestimados en 40% para el primer 
año y en 500% para el segundo; obte- 
niéndose un valor bajo de R2 en el que el 
modelo sólo explica un 42% de la varia- 
ción (Figura 7). Las sobreestimaciones 
que presenta el modelo, probablemente 
se deban a que los valores de los coefi- 
cientes de escurrimiento de los diferentes 
tipos de vegetación (curva numérica) no 
estaban del todo adecuadas a las condi- 
ciones de campo de la cuenca, así como 
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al modelo de evapotranspiración utilizado 
(Penman-Monteith). 

Para realizar una calibración satis- 
factoria se reajustaron varios parámetros 
del modelo, en especial las variables que 
tienen efecto en la evapotranspiración, la 
curva numérica y la capacidad de agua 
disponible en el suelo. Después de varios 
cambios y corridas del modelo, se llegó a 
una salida de datos que se ajustó mejor a 
la información de campo (R2 = 85%). El 
modelo se ajusta razonablemente a los 
datos en el primer año, sobre todo porque 
predice bien los picos máximos en los 
meses de mayor producción de agua. En 
el segundo año se observa una ligera 
sobreestimación pero la tendencia 
general de los meses de producción de 
agua se correlaciona satisfactoriamente 
(Figura 8). 



Mes (1994-1995) 



Figura 7. Salida y ajuste del modelo swat con las variables iniciales de la cuenca 

El Carrizal para los años 1994 y 1995. 
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Figura 8. Salida del modelo swat con el modelo calibrado para la cuenca El Carrizal 

con datos de los años 1994 y 1995. 


Fernández (1999) encontró en la 
cuenca agrícola El Tejocote en el Estado 
de México una buena eficiencia del 
modelo swat para predecir la producción 
de agua mensualmente con coeficientes 
de determinación (R2) de 0,80, 0,88 y 
0,61 en los tres años de medición. Por su 
parte, Benavides-Solorio (1998) obtuvo 
un valor de R2 de 0.91 en la cuenca El 
Carrizal en Jalisco. En relación a la 
predicción de la producción de sedi- 
mentos se encontró un valor de R2 
cercano a 0,65 en promedio en la cuenca 
El Tejocote y 0,78 en la cuenca El 
Carrizal. Además, Torres-Benites y cola- 
boradores (2004) encontraron valores de 
R2 de 0,94 y 0,84 para la producción 
anual de agua en dos años consecutivos 
en la cuenca El Tejocote. 


Producción de agua 

La tabla 2 muestra los resultados, 
por subcuenca, derivados del modelo 
calibrado. La producción de agua que 
genera el modelo se obtiene como una 
sumatoria de los escurrimientos superfi- 
ciales, escurrimientos subsuperficiales y 
escurrimientos laterales, a los cuales se 
les restan las pérdidas por transmisión. 

Producción de sedimentos 

La producción de sedimentos, como 
muestra el modelo, en las subcuencas se 
presentó durante la época de lluvias, prin- 
cipalmente en el periodo de junio a octubre 
en ambos años (Figura 9), comporta- 
miento similar al que presentaron los datos 
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Tabla 2. Componentes del balance hídrico por subcuenca de la cuenca Tapalpa. 


Subcuenca 

Precipitación 

(mm) 

1994 1995 

EV1 

(%) 

1994 

1995 

Producción 
de agua 
(%) 

ES2 Infiltración^ 
(%) (%) 

Tapalpa 

1 017 

812 

52 

78 

25 

12 

12 

El Zopilote 

1 017 

812 

54 

79 

22 

14 

8 

Las Ánimas 

1 017 

812 

60 

92 

9 

2 

3 

El Carrizal (Cal)4 

1 017 

812 

61 

95 

9 

2 

2 

El Carrizal (II) parte bajaS 

1 017 

812 

52 

73 

35 

30 

6 


1 Evapotranspiración. 

2 Se refiere al escurrimiento superficial y su aporte en la producción de agua. 

3 El agua que se infiltra del total del agua de lluvia. 

4 Esta subcuenca fue utilizada para calibrar el modelo, ya que en esta área se encuentran los instrumentos de 

medición. 

5 Esta subcuenca fue utilizada solo para fines comparativos, ya que la parte alta fue separada por la obra afora- 

dora. 


originales para la subcuenca El Carrizal 
(Benavides-Solorio, 1998). Los sedi- 
mentos son producto principalmente de la 
erosión causada por el impacto de las 
gotas de lluvia que ocasiona una disgrega- 
ción de los agregados del suelo y aumen- 
tada por el esfuerzo cortante del escurri- 
miento superficial. Las subcuencas que 
tuvieron la mayor producción de sedi- 
mentos fueron la parte baja de El Carrizal 
(II), Tapalpa y El Zopilote debido a que son 
aquellas que tienen mayor proporción de 
superficie agrícola. En contraste, la 
producción de sedimentos en las 
subcuencas Las Ánimas y El Carrizal (utili- 
zada para calibración) fue bastante baja, 
debido a la presencia de una alta cober- 
tura boscosa. La mayor producción de 
sedimentos se observó en el mes de 
octubre de 1994, periodo en el cual 
también se presentó la mayor producción 
de agua. En 1994, la mayor producción fue 
en octubre y en 1995 en septiembre. Los 
valores máximos fluctúan entre 10 y 17 


ton • ha-i -mes-l para las tres subcuencas 
antes descritas, mientras que para Las 
Ánimas y El Carrizal los valores máximos 
mensuales son muy bajos y del orden de 
0,1 ton -ha-i -mes-i, valores que se encuen- 
tran dentro de los límites de una cuenca 
forestal. 

En la Figura 9 se observa que la 
subcuenca baja de El Carrizal produce 
más sedimentos que las subcuencas 
Tapalpa o Zopilote, lo cual no necesaria- 
mente ocurre en la realidad por varias 
razones: (i) la subcuenca se delimitó arti- 
ficialmente; la parte baja no está sepa- 
rada de la parte alta; esto se hizo con el 
propósito de separar la parte instrumen- 
tada de la no instrumentada; (ii) un 
volumen considerable de agua proviene 
de la parte alta forestal (subcuenca 
Carrizal), por lo que los sedimentos se 
diluyen a mayor volumen de agua; (iii) 
debido a la proporción similar que guarda 
la cuenca completa El Carrizal con la 
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Mes (1994-1995) 


Figura 9. Producción de sedimentos por mes y subcuenca en la cuenca Tapalpa. 


subcuenca Las Ánimas, se estimaría que 
los valores de sedimento deberían tener 
valores similares en ambas subcuencas. 
Por lo tanto, los valores de sedimento de 
esta subcuenca deben manejarse con las 
debidas restricciones y como un antece- 
dente de lo que sucedería si se hicieran 
cambios en el uso del suelo, de uso 
forestal a uso agrícola o ganadero. 

Análisis de la condición hidrológica de 
la cuenca Tapalpa 

Para evaluar el aspecto hídrico y la 
funcionalidad de una cuenca, ésta debe 
estudiarse desde una perspectiva inte- 
gral. Para ello se toman en cuenta las 
características físicas de la cuenca, como 
son el suelo, la pendiente y la geología; 
las características de la vegetación 
(nativa e inducida), y áreas desnudas; los 
escurrimientos superficiales y subsuperfi- 
ciales, la evapotranspiración, el estado de 
las áreas ribereñas, y las fuentes y 
niveles de contaminación, entre otras. La 
cuenca Tapalpa tiene en sus partes altas 


un buen grado de conservación, debido al 
gran porcentaje de superficie cubierto por 
vegetación nativa. Los resultados del 
estudio permiten observar que las 
subcuencas con menores índices de 
producción de sedimentos son aquellas 
que tienen altos porcentajes de vegeta- 
ción forestal, como es el caso de las 
subcuencas El Carrizal y Las Ánimas. 

La precipitación de 914 mm en 
promedio de los dos años estudiados es 
una cantidad adecuada para la produc- 
ción de cultivos de clima templado-frío 
con buenos rendimientos. La producción 
anual de agua de 205 mm (22,5% del 
total de la precipitación) está balanceada 
por mitad en escurrimiento superficial y 
escurrimiento subsuperficial (Tabla 4), lo 
cual indica que 108 mm del total de la 
lluvia anual no se infiltran en el suelo 
(11,9% del total de la precipitación). Este 
escurrimiento superficial tiene la capa- 
cidad de arrastrar sedimentos y contami- 
nación a los arroyos y a la presa El Nogal; 
durante los periodos de tormentas, la 
producción de sedimentos se eleva en 
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magnitud por arriba de otros días en la 
temporada de lluvias. 

La vegetación de bosque combinada 
con suelos profundos proporciona una 
gran capacidad de almacenamiento de 
agua en las subcuencas, ya que aún 
durante los periodos de grandes 
tormentas los sedimentos y la contamina- 
ción son reducidos. En contraste, las 
subcuencas que presentan grandes 
áreas de agricultura o pastizal responden 
más rápido a las tormentas, incluyendo 
valores altos de escurrimientos superfi- 
ciales y grandes cantidades de sedi- 
mentos. Estas condiciones fueron detec- 
tadas por el modelo swat y siguen el 
mismo patrón con respecto a datos obte- 
nidos en campo en las subcuencas 
Tapalpa y El Zopilote (Benavides-Solorio 
etal., 2005a; 2005b). 

En la cuenca Tapalpa se presentan 
lluvias erosivas con alto potencial para 
producir sedimentos, por lo cual se debe 
tener una cubierta protectora. La intensidad 
máxima que se obtuvo en dos años de toma 
de datos fue de 37 mm en 30 minutos, la 
cual es una lluvia con alta capacidad erosiva 
(Benavides-Solorio, 1998), pero se tienen 
antecedentes de que pueden presentarse 
lluvias superiores a ésta. La manera de 
amortiguar estas lluvias es a través de la 
cubierta vegetal, sobre todo a nivel del suelo, 
lo que protegerá al mismo del impacto de las 
gotas de lluvia. En la zona no se han presen- 
tado niveles altos de deforestación, ya que 
los bosques están regulados por planes de 
manejo estrictos, pero en la última década 
se han presentado cambios de uso del suelo 
con fines recreativos, por lo que sería conve- 
niente realizar estudios sobre la dinámica 
del cambio de uso del suelo en la zona. 

La cubierta que reduce el impacto de 
las gotas de lluvia también retarda el escu- 
rrimiento superficial y permite que el agua 
se mantenga más tiempo sobre la super- 
ficie creando más oportunidades para que 
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ocurra la infiltración. Casi la mitad de la 
producción de agua en la cuenca Tapalpa 
se presenta como escurrimiento subsuper- 
ficíal debido a que la precipitación primero 
se infiltra en el suelo. La cubierta vegetal 
permite una infiltración inicial en el suelo y 
una liberación posterior del agua; el tiempo 
de retardo puede ser de horas o días 
dependiendo del contenido de humedad 
del suelo (Benavides-Solorio etal., 2005b). 
Si consideramos que la evapotranspira- 
ción en la cuenca es de 598 mm o un 
65.4% de la precipitación, es claro que 
existe la posibilidad de incrementar la 
producción de agua en la cuenca mediante 
el manejo adecuado de la cobertura 
vegetal (corta de árboles). Desafortunada- 
mente no se cuenta con una respuesta 
sobre cuál sería el método ni cuánto se 
debe cortar para incrementar la producción 
de agua, debido a que en México no se 
cuenta con antecedentes sobre el parti- 
cular. Es necesario realizar trabajos de 
investigación sobre cortas a diferentes 
intensidades y observar el efecto en la 
producción de agua. Aunque debe tenerse 
extremo cuidado para no aumentar la 
producción de sedimentos en los 
caudales. Lo anterior, es materia prima 
para futuras investigaciones. 

Finalmente, el volumen de agua 
disponible generado por la cuenca, se 
calculó mediante la lámina de riego de 
producción de agua, que es de 204,59 
mm (Tabla 3) y la superficie total de la 
cuenca (21,000 ha). De esta forma, el 
volumen de la oferta total asciende a 
42’963 900.0 m^-año-T Es necesario 
aclarar que los valores obtenidos para 
este trabajo fueron para un año promedio 
de lluvias y el valor obtenido puede 
aumentar o disminuir anualmente de 
acuerdo con la cantidad de precipitación 
anual. Para futuros estudios, si se desea 
conocer el comportamiento hidrológico 
para diferentes escalas, es necesario 
establecer parcelas experimentales de 
diferentes tamaños. 
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Tabla 3. Distribución del ciclo hidrológico anual en la cuenca Tapalpa. 


Variable 

Cantidad (lámina) 

Unidades 

Precipitación 

914,30 

mm 

Escurrimiento superficial (Q)"' 

107,98 

mm 

Escurrimiento lateraH 

18,37 

mm 

Escurrimiento subsuperficiaP 

78,98 

mm 

Evaporación 

5,24 

mm 

Recarga de acuíferos 

4,04 

mm 

Recarga total al suelo 

87,39 

mm 

Producción de agua 

204,59 

mm 

Percolación profunda 

89,18 

mm 

Evapotranspiración 

598,00 

mm 

Evapotranspiración potencial 

1320,00 

mm 

Pérdidas por transmisión 

0,73 

mm 

Producción de sedimentos 

19,60 

ton/ha 


1 Representan en conjunto la producción de agua. 


Implicaciones de uso del modelo swat 

En términos generales, y a pesar de 
su buena eficiencia, los modelos hidroló- 
gicos como el SWAT deben usarse sólo 
como apoyo para la toma de decisiones 
ya que, cuando se modelan escenarios 
futuros, pueden ocurrir modificaciones en 
campo que no fueron necesariamente 
consideradas en el o los modelos. 

Para tener mayor confiabilidad en la 
información que proporcionan los 
modelos, es necesario realizar calibra- 
ciones a los mismos con datos de campo, 
por lo cual se requiere la instrumentación 
de cuencas forestales para apoyar adicio- 
nalmente las decisiones de pagos por 
servicios ambientales hidrológicos. 

Por lo anterior, es necesario 
promover el establecimiento de cuencas 


experimentales en áreas forestales para 
obtener información confiable durante 
varios años. En este sentido, las cuencas 
pareadas han sido importantes como un 
método de campo eficiente para generar 
información con validez científica y en el 
corto plazo sobre la hidrología de una 
cuenca sometida a diferentes condiciones 
de manejo. Se trata de un método sencillo 
en su diseño pero que requiere mucho 
trabajo, en el cual se necesita medir prin- 
cipalmente el escurrimiento superficial en 
diferentes escalas de tiempo y espacio. 
Esta variable se mide durante dos etapas: 
calibración y post-tratamiento en las dos 
cuencas que se encuentran en condi- 
ciones ambientales similares, lo cual 
permite distinguir el efecto del tratamiento 
con respecto al efecto combinado de 
otros factores a lo largo del periodo de 
estudio, ya que una de las cuencas 
funciona como control del estudio. De 
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esta manera, a través de las mediciones 
en ambas cuencas es posible identificar 
el cambio relativo en la producción de 
agua con respecto a ios tratamientos apli- 
cados a la masa forestal y/o al suelo de la 
cuenca. 

Una de las primeras etapas en la 
implementación de un mercado de sah es 
realizar un diagnóstico de las caracterís- 
ticas socioeconómicas y biofísicas de la 
cuenca (Pagiola et al., 2003). Dentro de 
estas últimas, un aspecto importante es la 
estimación del balance hídrico de la 
cuenca. En este contexto, esta investiga- 
ción se puede utilizar para la aplicación 
de SAH con el uso del modelo swat, ya 
que permite conocer, además de la 
cantidad de agua superficial, la recarga 
de agua al subsuelo, la evapotranspira- 
ción y la producción de sedimentos, entre 
otros. Estos aspectos se consideran 
importantes para la toma de mejores 
decisiones para el pago de sah diferen- 
ciado. Terrenos forestales donde se capta 
agua, producto de los escurrimientos 
subsuperficiales, podrían tener un pago 
más alto, que aquellos terrenos donde 
existe mayor escurrimiento superficial, 
menor infiltración en el suelo y mayor 
arrastre de sedimentos. Una cuenca con 
alta capacidad de retención de agua 
muestra que la vegetación y los suelos se 
encuentran en condiciones estables y 
funcionales. La cobertura del suelo 
también es un indicador importante del 
comportamiento hidrológico que se puede 
esperar en una cuenca. Mientras más 
cubierto se encuentre el suelo se tendrá 
mayor protección, mayor infiltración y 
menor producción de sedimentos. 

Una gran ventaja al utilizar el modelo 
SWAT es que permite integrar diferentes 
suelos de las subcuencas, usos del suelo 
y tipos de vegetación, los cuales tienen 
respuestas diferentes en los escurri- 
mientos. Debido a que las cuencas de 
México pueden tener diferentes usos del 
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suelo, el modelo se ajusta a las varia- 
ciones del terreno y lo refleja en las 
respuestas. 

Una desventaja del modelo swat es 
que depende en gran medida de los 
valores de la curva numérica para el 
cálculo del escurrimiento superficial. A su 
vez, la curva numérica depende de la 
definición de la cobertura vegetal y el tipo 
de suelo, por lo tanto los límites deben 
estar definidos de manera precisa y 
actualizada en los mapas. Adicional- 
mente, el modelo swat no incluye en sus 
análisis a los caminos forestales, los 
cuales pueden ser productores de escu- 
rrimiento superficial y erosión en las 
zonas forestales. 

Por otra parte, es importante 
impulsar el desarrollo, adecuación y cali- 
bración de modelos existentes que 
permitan predecir la producción de agua 
en una cuenca, subcuenca o micro- 
cuenca en función del manejo que se le 
de a la cobertura vegetal; swat es uno de 
los modelos utilizados para ello; sin 
embargo, existen otros modelos como el 
AGWA (Automated Geospatial Watershed 
Assessment) (Semmens et al., 2004), el 
modelo hec-hms y varios otros que deben 
ser explorados a detalle para definir su 
utilidad para las cuencas forestales del 
país. Para el desarrollo de modelos 
propios es necesario contar con sufi- 
cientes datos que provengan de cuencas 
y parcelas experimentales, los cuales 
permitirían tener mayor seguridad para su 
calibración así como mayor confianza en 
su uso para las condiciones de los 
bosques de México. 

CONCLUSIONES 

El modelo swat permitió ajustar de 
manera satisfactoria los datos de la 
subcuenca El Carrizal, lo cual ayudó a 
conocer información no medida como la 
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recarga de mantos freáticos y evapotrans- 
pi ración. Las salidas del modelo swat en la 
producción de agua fueron similares a los 
datos obtenidos en campo, por lo cual 
puede utilizarse de manera confiable para 
conocer el ciclo hidrológico. 

La extrapolación de la calibración 
del modelo en la subcuenca El Carrizal 
hacia la cuenca Tapalpa proporcionó 
información adecuada en las salidas de 
agua y sedimento con respecto a los usos 
del suelo de cada subcuenca y ayudó a 
comparar la información para cada una 
de las cuencas, así como a proporcionar 
valores totales de la cuenca Tapalpa. 

Con base en lo anterior, las 
subcuencas que produjeron mayor 
cantidad de sedimentos fueron aquellas 
con mayor escurrimiento superficial y que 
también presentaban vegetación agricola 
y de pastizales: Tapalpa y El Zopilote, 
además de la parte baja de El Carrizal. 
Estas mismas subcuencas presentaron 
deterioro en la calidad de agua con base 
en las observaciones de Benavides-Solorio 
(2007) realizadas en el campo. Final- 
mente, debido al potencial turístico que 
impera en la zona, sería conveniente 
conservar y aumentar la condición arbo- 
lada de la cuenca Tapalpa. 
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Políticas que afectan la competitividad de la 
producción de madera aserrada en México 

Policies that affect the competitiveness of sawn wood 

production in México 

Martha E. Fuentes-Lópezi, J. Alberto García-Salazar2, Francisco 
J. Zannudio-Sánchez3 y J. Arturo Matus-Gardea2 

RESUMEN 

Un problema que enfrenta la producción nacional de madera aserrada es la pérdida de compe- 
titividad frente a la madera importada. Para el periodo 1993-1995, la producción nacional abastecía el 
80% del consumo nacional y para el 2001-2003 sólo el 3%. Con el objetivo de determinar los factores 
que podrían revertir esta tendencia, se usó un modelo de programación no lineal. Los resultados 
indican que una disminución de 30% en los costos de producción vía subsidios, reduciría las importa- 
ciones en 1,1%, estos es, 50,139 m^ aserrados, respecto al nivel observado en el periodo 2001-2003, 
y por otra parte incrementaría la producción en 2,1%. Los efectos de una reducción de 20% en los 
costos del transporte de las rutas internas nacionales causarían un aumento en la producción de 
3,8%, esto es, 90,500 m3 aserrados y un decremento en las importaciones de 2,2%, lo que sería igual 
a 102,034 m3 aserrados. Aplicar un arancel de 6 % sobre el precio internacional aumentará la produc- 
ción en 5,6%,, es decir, 132,355 m3 aserrados y reducirá las importaciones en 9,5%, que ascendería 
a 439,280 m3 aserrados, respecto a los niveles observados en 2001-2003. 

PALABRAS CLAVE: 

Demanda, importaciones, mercado, modelo de equilibrio espacial, oferta. 


ABSTRACT 


Mexican sawn wood is facing a lack of competitiveness with respect to imported wood. Over the 
period 1993-1995 average national production was able to supply 80% of the national consumption 
whereas by the period 2001-2003 oniy was capable of supplying 24%. In order to determine the factors 
that could reverl this tendency, a non linear programming model was used. Results indicate that if there 
were subsidies a 30% reduction in production costs would decrease imported wood in about 1,1%, 
which means 50,139 m3 of sawn wood, and it would increase the national production in 2,1%. A 20% 
reduction in transportation cost would lead to a rise in the national sawn wood production of 3,8%, this 
is 90 500 m3 in terms of sawn wood; consequently, imported sawn wood would decrease in 2,2%, 
which means 102 034 m3 of sawn wood. If a 6% tariff above the International price were to be imple- 
mented, the national production would increase in about 5,6%, this means 132 355 m3 of sawn wood. 
This tariff will help to decrease imported wood in about 9,5%, this is 306 923 m3 of sawn wood with 
respect to those of the 2001-2003 period. 

KEYWORDS: 

Demand, imports, market, spatial equilibrium model, supply. 
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INTRODUCCIÓN 

México ocupa el octavo lugar mundial 
y segundo en América Latina en extensión 
forestal. El 75% del territorio nacional 
corresponde a terrenos forestales, equiva- 
lente a 141,7 millones de ha, de los cuales, 
56,8 millones de ha están cubiertas por 
bosques y selvas. De esta superfide, 22 
millones tienen potencial de producción 
maderable sustentable con un volumen de 
38 millones de metros cúbicos rollo (Mm3r) 
(FAO, 1997; SEMARNAT, 2001 a 2003; Torres, 
2003). 

Los productos forestales pierden 
competitividad frente a la apertura comer- 
cial. En el caso de la madera, el aumento de 
las importaciones pone de manifiesto su 
falta de competitividad. El coeficiente de 
competitividad (razón de la producción 
sobre el consumo nacional de madera en 
rollo, p/c) ha descendido desde 1994, año 
de entrada en vigor del Tratado de Libre 
Comercio de América del Norte (tlcan). En 
el periodo promedio 1993-1995 la razón p/c 
fue 80% y para 2001-2003 disminuyó a 
34%. Como consecuencia, la razón de 
importaciones sobre consumo (m/C) fue 
poco más de 60% en los años 2002 y 2003. 

Entre las causas de la pérdida de 
competitividad destaca la liberalización 
comercial que ha afectado la producción, 
el consumo y las importaciones de 
madera. Debido a la disminución en los 
precios internacionales y nacionales, la 
producción de madera en rollo se redujo 
1,7 % anual, de 5.12 Mm3r en 1993, a 
4.32 Mm3r en 2003. La producción tuvo 
un decremento moderado, pero el 
consumo aumentó 10,7 % anual de 1993 
a 2003 (6.48 a 17.96 Mm3r), mientras que 
las importaciones pasaron de 0,816 a 
14.97 Mm3r (semarnat, 2001 a 2003). 

Otras causas que explican la 
pérdida de competitividad están relacio- 
nadas con la estructura de la industria 


forestal, la cual en su mayoría está 
formada por unidades de producción 
de pequeña escala con bajos niveles de 
producción, bajos rendimientos y escasa 
eficiencia de conversión de materia prima 
(Juárez y Hernández, 2002). La competi- 
tividad se puede analizar mediante la 
cadena productiva, la cual será competi- 
tiva si puede mantener ventajas compara- 
tivas que le permitan alcanzar, sostener y 
mejorar una posición en el mercado 
(SAGARPA, 2003). Por tanto, generalmente 
una cadena será competitiva si lo es en 
cada eslabón y una estrategia competitiva 
deberá contemplar los costos de produc- 
ción en cada eslabón. 

El bajo precio internacional de la 
madera aserrada resta competitividad a 
la producción nacional y estos productores 
tienen menores posibilidades de obtener 
ganancias. Para competir con las impor- 
taciones, la producción nacional de 
madera aserrada puesta en las regiones 
consumidoras debe venderse a un precio 
igual, o menor al de su competidor. Una 
disminución en ios costos de producción 
puede mejorar la competitividad, pues la 
madera nacional podría ser más barata a 
la importada. La disminución en el costo 
de producción podría ser consecuencia 
de una disminución en el precio de los 
insumos o por un incremento en la 
productividad. El menor costo de produc- 
ción permitiría vender a precios bajos 
manteniendo ganancias para el productor, 
y con ello, asegurar su permanencia en el 
mercado regional. 

Las zonas de mayor producción de 
madera se encuentran alejadas de las 
regiones de consumo. La mayor produc- 
ción de madera aserrada se concentra en 
el Norte, Occidente y Sureste del país, y 
la mayor demanda en el Centro. La movi- 
lidad del producto requiere el uso de 
transporte y cuando las distancias son 
mayores los costos de comercialización 
se elevan reduciendo las ganancias. La 
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topografía accidentada y la falta de una 
red caminera adecuada son factores que 
restan competitividad a la industria. Una 
disminución en los costos de transporte 
de las rutas comerciales internas mejo- 
raría la posición del producto nacional 
frente a las importaciones. 

Algunas políticas comerciales también 
permitirían mejorar la competitividad de 
la producción de madera. El estableci- 
miento de un arancel a las importaciones 
aumentaría el precio de la madera prove- 
niente del exterior haciendo más caro el 
producto extranjero. El precio de las 
materias primas e insumos usados en la 
producción, el grado de pulverización de 
la industria forestal, la tecnología, las 
condiciones obsoletas de la maquinaria y 
equipo, las tasas de interés y el régimen 
de propiedad ejidal que prevalece en 
México, son otros factores que tienen 
impacto sobre la competitividad. 

Por efecto de la liberación comercial 
se espera que las importaciones conti- 
núen incrementándose. La falta de 
competitividad de la producción nacional 
de madera hace indispensable proteger la 
actividad por las siguientes razones: a) la 
industria de la madera tiene encadena- 
mientos hacia atrás porque consume 
insumos que dejarían de producirse si la 
actividad desapareciera: b) la economía 
de miles de comunidades y millones de 
productores dependen de los ingresos 
generados en la actividad; c) el aumento 
de las importaciones y déficit comercial 
del sector silvícola se ha triplicado en los 
últimos años y; d) la industria del aserrío, 
segunda más importante dentro del 
sector forestal, genera empleos en 2058 
plantas (semarnat, 2003). 

OBJETIVOS 

El objetivo principal de la investiga- 
ción fue medir en qué magnitud respon- 


dería el consumo, la producción y las 
importaciones de madera aserrada a 
cambios en los costos de transporte, 
costos de producción y establecimiento 
de aranceles, y determinar cómo mejo- 
raría la competitividad del producto ante 
cambios en los factores antes mencio- 
nados. 

METODOLOGÍA 

El modelo 

Para alcanzar el objetivo se formuló 
y resolvió un modelo de equilibrio espa- 
cial de la madera aserrada. El modelo 
supone que existen regiones que comer- 
cian madera y cada región constituye un 
solo y distinto mercado. Las regiones 
están separadas por unidad física, pero 
no aisladas por los costos de transporte, 
y éstos son independientes del volumen, 
lo cual implica la inexistencia de econo- 
mías de escala. Conociendo las funciones 
de oferta y demanda del producto y los 
costos de transporte, el modelo permitió 
determinar la producción, el consumo, las 
importaciones y las movilizaciones de 
madera. 

Las razones que justifican el uso de 
un modelo de equilibrio espacial aplicado 
al mercado de la madera son: a) la 
producción y el consumo ocurren en lugares 
separados; b) productores, consumidores 
e importadores de madera responden 
de diferente manera a los estímulos de 
mercado y; c) la producción nacional no 
es suficiente para abastecer el consumo 
nacional aparente de madera, ya que 
más de 70% es abastecido con madera 
proveniente del exterior. 

Basados en Takayama y Judge 
(1971), la función objetivo del modelo de 
equilibrio espacial maximizó el Valor 
Social Neto (vsn), el cual es igual al área 
bajo la curva de demanda, menos el área 
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bajo la curva de oferta, menos el valor de 
las importaciones y menos los costos de 
transporte. Asumiendo s (s= 1,2...1Q4) 
regiones productoras de madera 
aserrada, d {d= 1 , 2 ... 85) regiones consu- 
midoras y m {m= 1,2... 96) puertos y fron- 
teras de internación de las importaciones, 
el modelo de programación cuadrática se 
puede expresar de la siguiente manera: 

10 9 

í=l m=l 
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md'^md . 


m=\ 


( 1 ) 


sujeta a: 
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4 Chihuahua, Durango, Michoacán, Jalisco, 
Oaxaca, México, Puebla, Guerrero, Veracruz e 
Hidalgo. 

5 Noroeste (integrada por Baja California, Baja 
California Sur, Sonora, Sinaloay Nayarit); Norte 
(Chihuahua, Coahuila, Durango, Zacatecas y 
San Luis Potosí); Noreste (Nuevo León y Tama u- 
lipas): Occidente (Jalisco, Michoacán, Colima, 
Aguascalientes y Guanajuato); Centro (Estado 
de México, Morolos, Hidalgo, Puebla, Tlaxcala y 
Querétaro); Sur (Oaxaca, Chiapas y Guerrero); 
Golfo (Veracruz y Tabasco) y; Península 
(Yucatán, Campeche y Quintana Roo). 

6 San Luis Río Colorado, Mexicali, Ciudad 
Juárez, Piedras Negras, Manzanillo, Mazatlán, 
Tuxpan, Tampico y Progreso. 




( 5 ) 


donde: Xq = intercepto de la función de 
demanda en la región d; yq = cantidad 
consumida de madera aserrada en la 
región d; coq = pendiente de la función de 
demanda de madera en la región d; 
ds = intercepto de la función inversa de 
oferta en la región s; ig = intercepto de la 
función de oferta en la región s; Xs = can- 
tidad producida de madera en la región s; 
í 7 s = pendiente de la función inversa de 
oferta de madera en la región s; = precio 
internacional de madera importada por 
el puerto rn; Xm = cantidad importada de 
madera aserrada por la frontera m; 
Psq = costo de transporte de madera de la 
región s a la región d; Xgq = cantidad de 
madera enviada de la región sala región d; 
Pmd = costo transporte de madera de la fron- 
tera m a la región d, x^q = cantidad de 
madera enviada de frontera m a región d. 

La ecuación (1) maximiza el vsn. La 
ecuación (2) establece que para cada 
región d, la cantidad consumida de 
madera aserrada (/d) es menor o igual, a 
la cantidad de producto que recibe la 
región d, proveniente de todas las regiones 
oferentes, tanto de zonas productoras del 
país, como de fronteras de internación. La 
ecuación (3) indica que la producción de 
madera aserrada en cada región produc- 
tora s (Xg) deberá ser mayor, o igual, al 
total de envíos de madera aserrada de las 
zonas productoras, hacia las regiones 
demandantes (x^d)- La ecuación (4) esta- 
blece que para cada frontera de interna- 
ción m, las importaciones de madera 
aserrada (Xm) deberán ser mayores o 
iguales a los envíos desde las fronteras de 
entrada de las importaciones a las distintas 
regiones de consumo (Xmd)- 

El modelo se validó comparando la 
producción, el consumo y las importa- 
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dones de madera aserrada observadas 
en el año promedio 2001-2003, con los 
valores arrojados por el modelo base. 
Una vez obtenido el modelo base se 
realizaron tres escenarios: a) reducción 
de 30% en el costo de producción de la 
madera aserrada; b) reducción de 20% 
en los costos de transporte de las rutas 
comerciales internas y, c) imposición de 
un arancel de 6% sobre las importa- 
ciones. 

Datos 

Para obtener la solución del modelo 
base se obtuvieron funciones de oferta y 
demanda por región, el precio interna- 
cional en los puntos de entrada de impor- 
taciones y los costos de transporte de las 
zonas productoras y puntos de entrada, a 
los centros consumidores. Las funciones 
de oferta y demanda se calcularon usando 
la elasticidad precio, las cantidades produ- 
cidas y demandadas y los precios al 
productor y consumidor. Las elasticidades 
precio de la demanda y oferta calculadas a 
nivel de región, provinieron de Fuentes 
(2006). La producción de madera aserrada 
se obtuvo de semarnat (2001 a 2003). Los 
precios anuales promedio al productor 
fueron proporcionados por aserraderos 
ubicados en distintas regiones del país^. 


7 Unión de Comunidades Forestales IXETO, S. A 
de C. V., Cámara de la Industria Mueblera del 
Estado de Jalisco, Cámara Nacional de la 
Industria Forestal Guadalajara, Madererías 
Carpicentro, S. A de C. V„ México, Comunidad 
Sta. María Zenzontepec, Oaxaca, Macosa- 
Madererías, Oaxaca, Empresa Forestal del 
Ejido “El Balcón”, Tecpan de Galeana, Guerrero, 
Empresa Forestal de San Juan Nuevo Paranga- 
ricutiro, Michoacán, Semarnat-Certificación 
Fitosanitaria, México, Semarnat-Estadísticas 
Forestales, México, Industrias Perroni, S. A. de 
C. V., Unión de Ejidos de Producción, Explota- 
ción, Comercialización e Industrialización Agro- 
pecuaria y Forestal de la Sierra Norte de 
Puebla, Fábrica de Tarimas, Empaques y 
Madera Dimensionada de Puebla, Madererías 
de Durango, Madererías de Chihuahua, Made- 
rerías del Estado de México. 


En la producción de madera aserrada se 
obtienen productos de diferentes gruesos, 
anchos y calidades que originan varios 
precios. Los precios de venta por metro 
cúbico de la madera aserrada al menudeo 
están referidos a la calidad “mill-run”, y 
fueron proporcionados por proveedores de 
madera y establecimientos ubicados en 
distintas regiones del país. 

El consumo anual aparente de 
madera por región se obtuvo multipli- 
cando el consumo por vivienda, por el 
número de viviendas. La información 
provino de inegi (2005) y semarnat (2001 
a 2003). 

El precio internacional de la madera 
aserrada se obtuvo de bancomext (2005) 
y de la Administración de Evaluación de 
Procedimientos de Aduanas del país. 
Dado que bancomext reporta datos en 
dólares, éstos se transformaron a pesos 
mexicanos mediante la tasa de cambio 
nominal reportada por el Banco de Infor- 
mación Económica del inegi para el 
periodo 2001-2003. 

En el cálculo de los costos de trans- 
porte se consideraron las tarifas promedio 
establecidas por varias empresas de 
transporte nacional^. Estas empresas no 
muestran tarifas uniformes y presentan 
un margen de variación que oscila de 
$6,00 a $40,00 por kilómetro recorrido. 
No obstante, entre 70 y 80% de los 
productores, distribuidores y consumi- 
dores de madera aserrada, utilizan 
camiones propios para realizar el trans- 
porte, pues representa un ahorro de 60%. 

El costo de transporte se obtuvo a 
través de la fórmula ^ Ppk)/CC, 

donde Psc/ es el costo de transporte de la 


8 Transportes Unidos de México (TUM), Centro 
de Distribución y Almacenaje de Transportes 
(DAT), EASO, CENTAURO, suvi, Julián de 
Obregón, Castores y iii Guerras. 
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zona de producción sala región de consumo 
d ($ por m3 aserrado), Dsq la distancia 
promedio óesad (km), Pp ^ el precio por kiló- 
metro de s a c/ ($) y CC la capacidad de 
carga del tráiler (m3 aserrados). 

La resolución del modelo se obtuvo 
utilizando el procedimiento minos, dise- 
ñado para problemas de optimización con 
funciones objetivo no lineales, escrito en 
el lenguaje de programación gams 
(General Algébrale Modeling Systems). 

RESULTADOS Y DISCUSIÓN 

El modelo fue validado en consumo, 
producción e importaciones por puerto y 
frontera. En la tabla 1 se observan los 
resultados con variaciones menores a 
10%, valores considerados permisibles 
entre los valores observados y estimados 
por el modelo base. Se observa que la 
solución del modelo subestima la produc- 
ción en -3,2% y sobreestima el consumo y 
las importaciones en 3,0 y 5,0%, respecti- 
vamente. En la producción de madera 
aserrada las diferencias entre los datos 
observados y estimados por el modelo 
base oscilaron entre -8,7% (Jalisco) y 
1,5% (Oaxaca). Para el consumo, las 
discrepancias entre los datos observados 
y los valores arrojados por el modelo base 
oscilaron entre: -7,3 % (Golfo) y 7,7% 
(Occidente). La mayor variación entre los 
volúmenes observadas y los estimados se 
obtuvo en el Golfo con 8,8%, en tanto que 
la menor variación se obtuvo en la frontera 
norte con -3,7%. Una vez validado el 
modelo se procedió a realizar escenarios 
que permitieran determinar los factores 
que podrían mejoran la competitividad de 
la producción nacional de madera 
aserrada. Es importante mencionar que 
dichos escenarios consideran ciertos 
supuestos ideales que están lejos de la 
problemática real. Por ejemplo, no se 
considera la capacidad real de producción, 
lo cual implica que la oferta puede incre- 


mentarse indefinidamente, sin considerar 
las limitaciones naturales de los bosques. 
El modelo solo consideró la distribución 
espacial en un año, sin contemplar las 
características temporales y el almacena- 
miento de madera. El país se dividió en 10 
regiones productoras y ocho regiones 
consumidoras, sin considerar la totalidad 
de las regiones. Por lo tanto, la veracidad 
de los resultados y el uso de éstos 
depende de la veracidad de los supuestos 
realizados, así como de la calidad de la 
información utilizada para obtener el 
modelo base. 

La principal fortaleza de la industria 
de la madera es el bajo costo de la mano 
de obra, su disponibilidad, calidad y 
además su productividad. Sin embargo, 
factores como los costos de materiales, 
falta de personal capacitado en todos los 
niveles, poco control de calidad, orienta- 
ción de la producción para satisfacer el 
mercado interno, así como la existencia 
de tecnología obsoleta, han originado 
una industria poco competitiva (Torres, 
2003). 

Bajo esta perspectiva, se consideró 
como primer escenario la reducción de 
30% en el costo de producción de madera 
aserrada en el año promedio 2001-2003. 
La respuesta por el lado de la oferta fue un 
incremento en la producción de 2,1% equi- 
valente a 50 1 38 m3 aserrados. Por la 
disminución en el costo de producción, el 
consumo no se ve afectado; sin embargo, 
el cambio en la producción ocasiona que 
las importaciones disminuyan en 1,1%, 
respecto a las observadas en 2001-2003. 
A nivel regional, el aumento en la produc- 
ción ocurre en todas las regiones produc- 
toras. Los resultados anteriores confirman 
lo señalado por Muñoz (2001), acerca de 
que la ventaja en los costos de producción 
(que pudiera reflejarse con mejoras en los 
productos y procesos), derivaría directa- 
mente en un aumento de la competitividad 
en el mercado de la madera aserrada. 
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Tabla 1. Validación del modelo de madera aserrada, 2001-2003. (m3 aserrados). 


Región 

Datos observados 

Modelo base 

Cambio 

Cambio 


2001-2003 



% 


m3 aserrados 




Producción 



Chihuahua 

511 372 

509 000 

-2,372 

-0,5 

Durango 

479 232 

459 167 

-20,065 

-4,2 

Michoacán 

300 307 

277 391 

-22,916 

-7,6 

Jalisco 

109 814 

100 303 

-9,511 

-8,7 

Oaxaca 

197 364 

200 294 

2,930 

1,5 

México 

65 899 

60 956 

-4,943 

-7,5 

Puebla 

80 831 

80 842 

11 

0,0 

Guerrero 

104 242 

101 739 

-2,503 

-2,4 

Veracruz 

38 965 

38 120 

-845 

-2,2 

Hidalgo 

36 740 

35 944 

-796 

-2,2 

Otros 

505 233 

488 182 

-17,051 

-3,4 

Nacional 

2 430 000 

2 351 938 

-78,062 

-3,2 



Consumo 



Noroeste 

753 505 

811 111 

57,606 

7,6 

Norte 

751 147 

793 333 

42,186 

5,6 

Noreste 

278 003 

288 056 

10,053 

3,6 

Occidente 

1 262 261 

1 360 000 

97,739 

7,7 

Centro 

2 229 563 

2 210 000 

-19,563 

-0,9 

Sur 

709 879 

758 000 

48,121 

6,8 

Golfo 

532 821 

493 846 

-38,975 

-7,3 

Península 

239 543 

244 225 

4,682 

2,0 

Nacional 

6 756 721 

6 958 571 

201,850 

3,0 



Importaciones 



Golfo 

2 238 141 

2 435 523 

197,382 

8,8 

Pacífico 

1 307 452 

1 360 000 

52,548 

4,0 

Frontera Norte 

842 247 

811 111 

-31,136 

-3,7 

Total 

4 387 840 

4 606 634 

218,794 

5,0 
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Las distancias entre zonas de produc- 
ción y consumo elevan los costos de trans- 
porte, restando competitividad al productor 
nacional de madera aserrada. Los costos de 
transporte disminuyen la ganancia de los 
productores. Un segundo escenario evaluó 
el impacto que genera una reducción de 
20% en los costos de transporte. Por el lado 
de la oferta, la respuesta fue un incremento 
en la producción de 3,8% (90 500 m3 
aserrados), y por el lado de la demanda, se 
presentó un ligero decremento en el 
consumo de 0,2 % (-11 533 m3 aserrados). 
Las importaciones se redujeron en 2.2% 
(-102 034 m3 aserrados), respecto al nivel 
observado en 2001-2003. 

Las regiones más beneficiadas por la 
disminución en los costos de transporte 
serían: Chihuahua, Durango y Michoacán, 
en donde la producción se incrementaría 
en 30 000, 1 8 500 y 4 609 m3 aserrados, 
respectivamente, y solo el Norte y Sur se 
verían afectadas por el lado del consumo. 
Las importaciones realizadas por el Golfo 
de México se verían disminuidas. 

Una disminución en los costos de 
distribución y producción podría lograrse a 
través del aumento de la competitividad de 
los productos que se usan como insumos en 
el proceso de producción y comercialización 
de madera aserrada. Una disminución en 
los precios del diesel bajaría los costos de 
transporte con los efectos sobre la produc- 
ción y las importaciones que ya fueron anali- 
zados. Aunque actualmente se subsidia el 
precio del diesel en la fase primaria de 
producción, éste no ha sido suficiente para 
impactar sobre la cadena productiva. Por las 
ineficiencias que podrían generarse, no es 
recomendable un subsidio al distribuidor de 
madera aserrada, mejor se recomienda que 
el gobierno implemento las medidas nece- 
sarias para hacer competitiva la producción 
del diesel que se produce en el país. 

Otro forma en que podría disminuirse 
el costo de producción sería aumentando 


la productividad. Los volúmenes de 
producción de madera aserrada son muy 
bajos, principalmente por la falta de una 
infraestructura tecnológica y de innova- 
ción. El apoyo de actividades tecnoló- 
gicas, de innovación e investigación, con 
el propósito de incrementar la producti- 
vidad, disminuiría el costo medio de 
producción, mejorando la competitividad. 

Las empresas de la industria forestal 
enfrentan costos elevados de la madera 
aserrada, costos que se atribuyen a las 
operaciones de colecta, transporte y 
distancias de acarreo. Estos aspectos 
ejercen una influencia importante en la 
capacidad de las industrias para competir 
y ampliarse. Aunque la ventaja en los 
gastos de mano de obra y distribución al 
mercado están a favor de los produc- 
tores; el saldo del costo global es tal, que 
no podrían absorber los incrementos en 
los costos de transformación de la 
madera y mantenerse en una situación 
competitiva, a menos que se estable- 
cieran barreras arancelarias como una 
medida de protección (fao, 2005). 

Los efectos de un arancel sobre las 
compras externas elevan el precio 
doméstico de la madera importada e 
incrementan la producción nacional. El 
consumo doméstico se contraería origi- 
nando que las importaciones disminuyan 
también. El efecto del arancel dependerá 
de la sensibilidad de la oferta y la 
demanda de madera aserrada a cambios 
en el precio. Debido a la elasticidad 
precio que caracteriza a la oferta y la 
demanda de madera aserrada, existen 
elementos para suponer que un cambio 
en el precio de las importaciones tendrá 
efectos importantes sobre la cantidad 
importada, y como medida de política 
para proteger al productor doméstico y 
provocar un efecto significativo en el 
volumen de importación, se requerirá 
imponer una tasa arancelaria a las impor- 
taciones de madera aserrada. 
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Tabla 2. Factores que afectan la competitividad de la madera aserrada (m3 aserrados). 


Región 

Modelo 

base 

Disminución Disminución 
en 30% en 20% 

en CPU en CT§ 

Establecimiento 
de un arancel 
de 6 % 


Cambio 



m3 aserrados 






Producción 





Chihuahua 

509 000 

514 400 

539 000 

543 560 

1,1 

5,9 

6,8 

Durango 

459 167 

464 417 

477 667 

487 967 

1,1 

4,0 

6,3 

Michoacán 

277 391 

282 087 

282 000 

291 087 

1,7 

1,7 

4,9 

Jalisco 

100 303 

104 621 

102 970 

105 076 

4,3 

2,7 

4,8 

Oaxaca 

200 294 

204 882 

203 235 

209 575 

2,3 

1,5 

4,6 

México 

60 956 

65 147 

61 191 

63 272 

6,9 

0,4 

3,8 

Puebla 

80 842 

85 074 

81 284 

84 158 

5,2 

0,5 

4,1 

Guerrero 

101 739 

106 043 

103 246 

106 304 

4,2 

1,5 

4,5 

Veracruz 

38 120 

42 248 

38 488 

39 380 

10,8 

1,0 

3,3 

Hidalgo 

35 944 

40 066 

36 084 

37 045 

11,5 

0,4 

3,1 

Otros 

488 182 

493 091 

517 273 

516 869 

1,0 

6,0 

5,9 

Nacional 

2 351 938 

2 402 076 

2 442 438 

2 484 293 

2,1 

3,8 

5,6 




Consumo 





Noroeste 

811 111 

811 111 

811 111 

784 178 

0 

0,0 

-3,3 

Norte 

793 333 

793 333 

780 000 

754 933 

0 

-1,7 

-4,8 

Noreste 

288 056 

288 056 

288 056 

278 456 

0 

0,0 

-3,3 

Occidente 

1 360 000 

1 360 000 

1 360 000 

1 297 000 

0 

0,0 

-4,6 

Centro 

2 210 000 

2 210 000 

2 210 000 

2 105 000 

0 

0,0 

-4,8 

Sur 

758 000 

758 000 

759 800 

726 444 

0 

0,2 

-4,2 

Golfo 

493 846 

493 846 

493 846 

466 431 

0 

0,0 

-5,6 

Península 

244 225 

244 225 

244 225 

239 206 

0 

0,0 

-2,1 

Nacional 

6 958 571 

6 958 571 

6 947 038 

6 651 648 

0,0 

-0,2 

-4,4 



Importaciones 





Golfo 

2 435 523 

2 385 384 

2 333 489 

707 498 

-2,1 

-4,2 

-71,0 

Pacífico 

1 360 000 

1 360 000 

1 360 000 

2 675 678 

0,0 

0,0 

96,7 

Norte 

811 111 

811 111 

811 111 

784 178 

0,0 

0,0 

-3,3 

Total 

4 606 634 

4 556 495 

4 504 600 

4 167 354 

-1,1 

-2,2 

-9,5 


Costo de producción. 


§ Costos de transporte. 
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El tercer escenario consistió en 
aplicar un arancel de 6% sobre el precio 
de importación de la madera aserrada. 
Los resultados señalan que el estableci- 
miento del arancel causará un incremento 
en la producción en 5,6%, respecto de la 
producción registrada en 2001-2003. Esto 
significa que la producción aumentará en 
132 355 m3 aserrados. 

La medida proteccionista disminuirá el 
consumo en 4,4% respecto al nivel obser- 
vado en 2001-2003, equivalentes a 306 
923 m3 aserrados. Las importaciones 
disminuirían en 9,5%, respecto al nivel 
observado en el año de análisis; esto es, 
439 280 m3 aserrados. Las regiones bene- 
ficiadas por esta política serían las mayores 
oferentes como Chihuahua, Durango y 
Michoacán, que verían incrementada su 
producción en 34 560, 28 800 y 1 3 696 m3 
aserrados, respectivamente. Las regiones 
consumidoras más afectadas serían el 
Occidente y el Centro. 

Si se considera el actual escenario 
de apertura comercial y eliminación de 
barreras proteccionistas, la política 
de eliminación de aranceles no sería 
viable. Sin embargo, los fuertes efectos 
que la liberación comercial ha tenido 
sobre el mercado de la madera, y la 
importancia que el producto tiene para los 
demás sectores de la economía, justifican 
la propuesta de establecer un arancel 
para proteger la actividad. 

CONCLUSIONES 

La producción y el consumo de 
madera aserrada son sensibles al esta- 
blecimiento de aranceles, y en menor 
medida a cambios en los costos de 
producción y transporte. Un desplaza- 
miento importante de las importaciones 
podría lograrse a través de la instrumen- 
tación de políticas conjuntas que contem- 
plarán el establecimiento de una tarifa 


arancelaria mayor al 6% y la disminución 
en los costos asociados a la producción y 
distribución del producto. Debido a los 
encadenamientos que el sector forestal 
tiene hacia atrás, a las miles de comuni- 
dades que dependen de las actividad 
silvícola y a la fuerte presión que sobre la 
balanza de pagos ejercen las importa- 
ciones de madera, es recomendable la 
instrumentación de dichas políticas, pues 
permitirían mejorar la competitividad de la 
producción nacional del producto, aunado 
a un esfuerzo por crear y aumentar las 
ventajas competitivas que pudieran 
conducir hacia una real consolidación en 
los mercados internacionales. 
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lARTICULO DE INVESTIGACION 


Caracterización de medios de crecimiento 
compuestos por corteza de pino y aserrín 

Growth media characterization based on pine bark 

and sawdust 

Tito Sánchez-Córdova'', Arnulfo Aldrete2*, 

Víctor M. Cetina-Alcalá2 y Javier López-Uptori2 

RESUMEN 

En México, la producoión de plantas forestales en contenedores rígidos se realiza utilizando principal- 
mente como medio de crecimiento la turba, la vermiculita y la agrolita (6:3:1 , mezcla estándar). El aserrín y 
la corteza de pino son desechos de la industria forestal que pueden ser usados como sustratos. Este trabajo 
tuvo como objetivo caracterizar cinco mezclas de sustratos, compuestos de corteza, aserrín, turba, vermicu- 
lita y agrolita. La evaluación de las mezclas incluyó la caracterización física (densidad aparente, porosidad y 
curva de retención de humedad) y química (pH, conductividad eléctrica y capacidad de intercambio catió- 
nico). La densidad aparente de las mezclas compuestas de corteza y aserrín varió de 0,31 a 0,51 g cm-3 y 
la estándar fue de 0,15 g cm-3 La porosidad total de 77% a 91% y la estándar fue de 84%. La porosidad de 
aireación varió de 8% a 20% y la estándar fue de 8%. El agua fácilmente disponible varió de 4,6% a 9,9% y 
la estándar presentó un valor más elevado (17,7%). El pH de las mezclas compuestas de corteza y aserrín 
varió de 4,52 a 7,70 y la estándar fue de 3,87. La conductividad eléctrica de 1 ,84 a 2,69 ds m-i y la estándar 
1,13ds nri. La capacidad de intercambio catiónico de 43,6 a 89,4 meq/IOOgy la estándar 60,4 meq/100 g. 
La mezcla 4 (20% de corteza + 80% de aserrín) es semejante a la mezcla estándar por sus características 
físicas, pero no en sus propiedades químicas. Esta mezcla es un sustrato prometedor en la producción de 
plantas forestales. 

PALABRAS CLAVES: 

Aserrín, corteza, propiedades físicas, propiedades químicas, sustratos, turba. 


ABSTRACT 

Forest nursery production in containers in México mainly uses peat moss. vermiculite and agrolite 
(standard mixture). Sawdust and pine bark are subproducts of the forest industry which can be used as subs- 
traías. The objective of this work was the characterization of five mixtures using pine bai1<, sawdust, peat 
moss, vermiculite. Evaluation included physicai (bulk density, total porosity, air porosity and easily available 
water) and Chemical (pH, electrical oonductivity and catión exchange capacity) characterization. Bulk density 
variad from 0,31 to 0,51 g m-3 while the standard mixture was 0,15 g m-3. Total porosity varied from 77 to 
91% compared to the standard mixture with 84%. Ar porosity varied from 8 to 20% compared to the stan- 
dard mixture with 8%. Easily available water varied from 4,8 to 9,9% compared to the standard mixture with 
17,7%. pH of the mixtures with bark and sawdust varied from 4,52 to 7,70 while the standard mixture was 
3,87. Electrical conductivity varied from 1,84 to 2,69 ds m-l while the standard mixture was 1,13 ds m-i. 
Catión exchange capacity varied from 43,6 to 89,4 meq/100 g while the standard mixture was 60,4 meq/100 
g. Physicai properties of mixture 4 (20% of bark + 80% of sawdust) are similar to standard mixture, but not 
the Chemical properties. This mixture is a promising substrate for forest nursery production. 

KEY WORDS; 

Bark, Chemical properties, peat moss, physicai properties, sawdust, substraías. 
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INTRODUCCIÓN 

La calidad de las plantas que se 
producen en vivero depende, entre otros 
factores, de la adecuada selección de los 
sustratos para la preparación de medios 
de crecimiento. Una mezcla adecuada 
debe tener propiedades físicas y 
químicas que permitan la disponibilidad 
oportuna de los nutrimentos y el agua 
(Burés, 1997). El medio de crecimiento es 
uno de los factores que influye directa- 
mente en la calidad y costo de producción 
de las plantas en vivero, por ello se deben 
buscar opciones que reduzcan esos 
costos y garanticen la calidad de la planta 
(Arteaga etal., 2003). 

Actualmente en México se usa como 
sustrato principal, en la producción en 
contenedores rígidos, una mezcla de 
turba, agrolita y vermiculita (principalmente 
en proporciones de 60:30:10). Estos mate- 
riales son importados a altos costos, a 
pesar de que existen otros materiales que 
pueden ser útiles como sustratos y entre 
los que se encuentran el aserrín y la 
corteza de pino. El aserrín y la corteza son 
residuos del proceso de aserrío y pueden 
llegar a ser un problema en la industria de 
la madera, debido a que ocupan mucho 
espacio en la línea de producción después 
del aserrío (Juárez et al., 2001). Estos 
materiales son más baratos (hasta 70% 
menos) que la turba, vermiculita y agrolita, 
y también tienen características apro- 
piadas para reducir la actividad de hongos 
fitopatógenos y mejorar la porosidad 
(Landis et al., 1990). La corteza aumenta 
su capacidad de intercambio catiónico 
cuando está composteada. El aserrín de 
pino puede presentar algunos problemas 
de fitotoxicidad cuando se usa crudo, pero 
el problema se corrige con el lavado del 
mismo o con el proceso de descomposi- 
ción (Boodley, 1998). 

Cada mezcla utilizada en la produc- 
ción de planta forestal tiene propiedades 
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de densidad, porosidad, conductividad 
eléctrica, capacidad de intercambio catió- 
nico y pH que le proporcionan caracterís- 
ticas distintas de desarrollo y crecimiento 
a la planta (Burés, 1997). Algunos autores 
sugieren que el pH del medio de creci- 
miento de las especies coniferas debe 
estar entre 5,5 y 6,5, un pH por abajo o 
por arriba de estos rangos puede causar 
problemas en la germinación y en el 
crecimiento de la planta (Peñuelas y 
Ocaña, 1996). 

El presente trabajo tuvo como obje- 
tivo la caracterización de cuatro mezclas 
de sustrato compuestas de corteza y 
aserrín, con base en las propiedades 
físicas y químicas, y compararlas con la 
mezcla estándar utilizada en la produc- 
ción en contenedores, en muchos viveros 
forestales tecnificados de México. 

MATERIALES Y MÉTODOS 

La preparación de las mezclas se 
realizó en el vivero de Posgrado Forestal 
en las instalaciones del Colegio de 
Posgraduados en Montecillo, Estado de 
México. Éste se localiza en las coorde- 
nadas 19° 29’ LN y 98° 54’ LO, a una 
altitud de 2240 m (Reyes, 2005). Los 
materiales usados para las mezclas 
fueron: corteza de pino composteada, 
aserrín de pino, turba, agrolita y vermicu- 
lita (Tabla 1). La corteza utilizada fue de 
Pinus douglasiana y se obtuvo de la 
región sur del estado de Jalisco, donde 
es composteada y se vende en forma 
comercial., El aserrín fue de Pinus patula 
SchI. et Cham., recolectado crudo direc- 
tamente del aserradero del ejido Pueblo 
Nuevo, en el municipio de Zacatlán, 
Puebla, en donde se ha utilizado para la 
producción de especies forestales en 
vivero. La turba, la agrolita y la vermicu- 
lita son materiales importados de 
Canadá y se venden comercialmente en 
México. 
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Tabla 1 . Proporción de los componentes en las mezclas evaluadas. 


Mezcla 

Corteza 
de pino 
composteada 

Aserrín 

{%)(%) 

Turba 

(%) 

Agrolita 

{%) 

Vermiculita 

{%) 

1 

100 





2 

80 

20 




3 

60 

40 




4 

20 

80 




5* 



60 

30 

10 


* = Mezcla testigo 


El análisis de las propiedades físicas 
y químicas se realizó en 2006 en el Labo- 
ratorio de Física de Suelos del Colegio de 
Posgraduados. Para la determinación de 
las propiedades físicas de los sustratos 
se utilizaron cuatro muestras (repeti- 
ciones) de cada uno de ellos, tomadas al 
azar, con un volumen promedio de 
137 cm3, equivalente al volumen de un 
contenedor para la producción de planta. 
La densidad aparente se calculó con 
base en la fórmula Da = peso del sustrato 
seco (g) / volumen total (cm3). El peso del 
sustrato se obtuvo con una balanza analí- 
tica después de secarlo en estufa durante 
72 horas. La porosidad de retención de 
agua (prh), porosidad de aireación (pa) y 
porosidad total (pt), se calcularon con el 
procedimiento descrito por Landis et al. 
(1990). 

De una muestra de cinco litros de 
cada sustrato, se llenaron cuatro contene- 
dores con un volumen de 137 cmS; poste- 
riormente, a los contenedores se les 
agregó agua hasta saturarlos y luego se 
colocó un envase en la parte inferior para 
colectar el agua drenada. Con estos 


datos se calculó la porosidad. Las propie- 
dades químicas que se analizaron se 
basaron en el manual de laboratorio 
propuesto por Palmer y Troeh (1989) y 
van Reeuwijk (1999); se analizó el pH por 
el método del potenciómetro, la conducti- 
vidad eléctrica (ce) por el método del 
conductímetro y la capacidad de inter- 
cambio catiónico (cic) por el método de 
acetato de amonio. 

Con base en la curva de retención 
de humedad, se determinó la curva de 
liberación de agua para cada mezcla, 
usando el método de batería de embudos 
(De Boodt et al., 1974). La curva está 
compuesta de agua difícilmente dispo- 
nible (add), agua fácilmente disponible 
(afd), agua de retención (ar), capacidad 
de aire (ca), espacio poroso total (ept) y 
materia sólida (ms). 

Los datos de las propiedades físicas 
y químicas de las mezclas de sustrato 
(excepto los de la curva de retención de 
humedad) se analizaron mediante un 
análisis de varianza (anova) de dos vías 
de clasificación, más la interacción. 
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RESULTADOS Y DISCUSIÓN 
Propiedades químicas de las mezclas 

Existe variación en las propiedades 
químicas de las cinco mezclas, depen- 
diendo de los componentes en las 
mismas (Tabla 2). El pH de todas las 
mezclas está debajo de un valor de 5. La 
mezcla testigo fue la más ácida con un 
valor promedio de 3,87. En las mezclas 
que contienen corteza y aserrín, el pH es 
también ácido y varía de 4,25 en la 
mezcla 2 hasta 4,70 en la mezcla 4. 
Landis et al. (1990) reportan que con un 
pH menor a 4,0 se pueden manifestar 
problemas de enfermedades de raíz. Con 
valores de pH entre 5,5 a 6,5 el problema 
de fungosis se reduce y con pH mayores 
a la neutralidad pueden presentarse 
problemas por Fusaríum. En el medio de 
crecimiento puede haber variaciones en 
el pH al final de la producción, debido a 
las prácticas de cultivo (fertilización y 
riego), en la producción de plantas fores- 
tales puede haber un incremento de 0,5 a 
1,0 unidades de pH (Landis et al., 1990) 
por lo cual incluso la mezcla testigo 
podría alcanzar niveles adecuados de pH 
al momento de la producción. 
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La conductividad eléctrica (ce) de la 
mezcla testigo es la más baja con un valor 
promedio de 1 ,1 3 ds m-'' y ios valores más 
altos pertenecen a las mezclas 1 y 3 con 
2,69 y 2,68 ds m-’i, respectivamente. De 
acuerdo con Moreno (1970), no se 
detectan efectos detrimentales en el creci- 
miento y rendimientos de las plantas con 
una CE menor a 2,0 ds m-i, pero de 2,0 a 
4,0 ds m-i se restringen los rendimientos 
de algunas plantas sensibles. Bunt (1988) 
indica que en sustratos ricos en materia 
orgánica la ce debe ser de 0,75 a 1,99 
ds m-i y que de 2,00 a 3,49 ds m-i las 
plantas crecen satisfactoriamente, pero si 
son sensibles a ce por arriba de 2,0 ds m-i 
reducirán su crecimiento. También señala 
que en una ce de 3,5 a 5,0 ds m-i la sali- 
nidad es muy elevada para la mayoría de 
las plantas, y que sólo plantas vigorosas 
resisten esta condición. 

La capacidad de intercambio catió- 
nico (cic) de las mezclas 1 y 5 (testigo) 
son semejantes con 61,2 y 60,4 meq/ 
lOOg. Al respecto, Peñuelas y Ocaña 
(1996) indican que la corteza de pino 
tiene una cic de 52,6 meq/IOOg. Estos 
autores también mencionan que la cic 
aumenta cuando incrementa el pH y que 


Tabla 2. Propiedades químicas de las mezclas evaluadas en el laboratorio. 


Mezclas Porcentaje de componentes Propiedades químicas 



de las mezclas 

pH 

relación 1:2 

CE CIC 

(ds m-i) (meq/IOOg) 

1 

Corteza 100% 

4,69 

2,69 

61,2 

2 

Corteza 80% + Aserrín 20% 

4,52 

1,84 

89,4 

3 

Corteza 60% + Aserrín 40% 

4,63 

2,68 

87,6 

4 

Corteza 20% + Aserrín 80% 

4,70 

2,58 

43,6 

5* 

Turba 60% + Agrolita 30% + Vermiculita 10% 

3,87 

1,13 

60,4 


*= mezcla testigo, ce= conductividad eléctrica y cic= capacidad de intercambio catiónico. 
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ambos aumentan cuando incrementa el 
contenido de materia orgánica, además 
amortigua los efectos de pH y la sali- 
nidad. En el caso de las mezclas anali- 
zadas no se incrementó la cic cuando 
aumentó el pH. Landis et al. (1990) 
indican que la turba y la vermiculita tienen 
valores de cic con 180 meq/IOOg y 82 
meq/IOOg, respectivamente, mientras 
que la perlita (agrolita) tiene un valor 
promedio de sólo 3,5 meq/IOOg. La 
mezcla de 80% de corteza + 20% de 
aserrín y la mezcla de 60% de corteza + 
40% de aserrín tienen valores similares en 
cic. A mayores proporciones de corteza y 
menores de aserrín, aumentó la cic. La 
mezcla 4 con porcentaje de aserrín del 
80% tiene la cic más baja (43,6), lo que 
puede ser debido a que el aserrín crudo 
no tiene los suficientes microorganismos 
que descompongan la materia orgánica y 
por tanto su cic es menor. 

Landis et al. (1990) señalan que los 
materiales orgánicos no composteados, 
como el aserrín o la corteza, limitan la dispo- 
nibilidad de nitrógeno, ya que los orga- 
nismos que descomponen la materia orgá- 
nica utilizan este nutrimento, y por tanto su 


cic se reduce. Entre mayor sea la cic es 
mejor para una buena acumulación de 
cationes disponibles para las plantas, y si la 
mezcla tiene una cic reducida hay deficien- 
cias de nutrimentos (Landis et al., 1990). 


Propiedades físicas de las mezclas 

En cuanto a las características 
físicas de las mezclas, las densidades 
reales en seco (drs) oscilan entre 1 ,09 y 
1,51 g cm-3 donde la mezcla 1 (100%) de 
corteza es la más alta y la mezcla testigo 
es la de menor valor (Tabla 3). Las densi- 
dades reales en húmedo (drh) están 
entre 1,03 y 1,15 g cm-3. La mezcla 1 fue 
nuevamente la más alta y la mezcla 
testigo la más baja. Burés (1997) 
menciona que la densidad real para la 
corteza de pino triturada es de 1,84 g cm-3. 
Landis et al. (1990) señalan que las 
densidades de los minerales son de 2,6 g 
cm-3 y de la materia orgánica es de 1 ,55 
g cm-3. Sin embargo, las densidades de 
las mezclas pueden variar por el grado de 
compresibilidad y del arreglo de las partí- 
culas. La densidad disminuyó conforme 
se reduce la proporción de corteza. 


Tabla 3. Propiedades físicas de las mezclas evaluadas en el laboratorio. 


Mezcla 



Propiedades físicas 





DRS 

DRH 

DAS 

DAH 

PT 

PA 

PRH 


g cm-3 

g cm-3 

g cm-3 

g cm-3 

% 

% 

% 

1 

1,51 

1,15 

0,51 

1,04 

70 

8 

62 

2 

1,48 

1,14 

0,43 

1,01 

77 

9 

68 

3 

1,39 

1,14 

0,40 

1,01 

78 

9 

69 

4 

1,40 

1,10 

0,31 

0,94 

91 

20 

71 

5* 

1,09 

1,03 

0,15 

0,82 

84 

8 

76 


* = Mezcla testigo, drs= densidad real en seco, drh= densidad real en húmedo, das= densidad aparente en seco, 
DAH= densidad aparente en húmedo, pt= porosidad total, pa= porosidad de aireación y prh= porosidad de retención 
de humedad. 


Mezcla 1 (corteza al 100%), mezcla 2 (corteza 80% y aserrín 20%), mezcla 3 (corteza 60% y aserrín 40%), mezcla 
4 (corteza 20% y aserrín 80%) y mezcla 5 (turba 60%, agrolita 30% y vermiculita 10%). 
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Tensión en cm de columna de agua 


EPT 
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b) Corteza 80% + aserrín 20% 
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c) Corteza 60% + aserrín 40% 
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d) Corteza 20% + aserrín 80% 
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e) Turba 60% + agrolita 30% + vermiculita 10% 


Figura 1. Curvas de liberación de agua de las cinco mezclas analizadas en el estudio. 

ADD = agua difícilmente disponible, ar = agua de reserva, afd = agua fácilmente 
disponible, ca= capacidad de aireación, ms = materia sólida y ept = espacio poroso total. 
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También se aprecia que las densi- 
dades aparentes en seco (das) de las 
mezclas oscilan entre 0,15 y 0,51 g cm-3, 
donde la mezcla 1 es la más alta y la 
mezcla 5 la más baja (Tabla 3). En el caso 
de la densidad aparente en húmedo (dah) 
los valores están entre 0,82 y 1 ,04 g cm-3 
y corresponden a la mezcla 1 que es la 
más alta y la mezcla 5 la más baja (Tabla 
3). La DAH es importante para el manejo 
de los contenedores en vivero o para el 
traslado de las plantas al campo. Mezclas 
con bajas dah permiten manejar mayor 
número de contenedores porque tienen 
menor peso y llevar más carga de 
plantas. La das aumentó con respecto a 
la DAH en las mezclas con mayor propor- 
ción de corteza. Por otra parte, casi en 
todas las mezclas el valor de la dah es 
semejante, ya que la máxima diferencia 
entre el valor más alto y el más bajo es de 
0,2 g cm-3. ei incremento de la das a dah 
en todas las mezclas fue mayor de 0,5 g 
cm-3, es decir, aumentan por lo menos en 
100% de peso de seco a peso húmedo. 
Esta diferencia indica que habrá alta 
capacidad de retención de humedad en 
todas las mezclas. 

La porosidad total (pt) en la mezcla 
4 es la más alta con sólo 9% de material 
sólido, mientras que la mezcla 1 presenta 
la menor pt y 30% de material sólido 
(Tabla 3). La pt en las mezclas 2 y 3 son 
similares, con material sólido de 23% y 
22%, respectivamente. La pt en la 
mezcla testigo es de 84% y 16% de mate- 
rial sólido. Landis et al. (1990) señalan 
que muchos autores recomiendan de 60 
a 80% de porosidad total. De acuerdo con 
esto, tres de las mezclas probadas están 
dentro de esos márgenes. Con respecto a 
la porosidad de aireación (pa), todas las 
mezclas están entre 8% y 20% (Tabla 3); 
la mezcla 4 es la que presenta una mayor 
PA con relación a las otras. Con respecto 
a la porosidad de aireación, Havis y 
Hamilton (1976) recomiendan que la pa 
debe ser de 20% a 25%; Whitcomb 
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(1988) de 25% a 35% y Puustjarvi y 
Robertson (1975) de 45% a 50%. Los 
valores de pa presentes en la mayoría de 
las mezclas indican que existen pocos 
macroporos y que las raíces de las 
plantas tienen poca aireación. Lo que 
supone que el tamaño de las partículas 
es fino, hay microporos y la retención de 
humedad es alta. Lo que puede indicar 
posibles problemas de podrición de 
raíces por falta de oxígeno. Una razón de 
los bajos porcentajes de poros de airea- 
ción es la variabilidad de la granulometría 
de los sustratos usados, ya que las partí- 
culas menores de los sustratos entran a 
los macroporos y disminuye la pa. Lo 
anterior se corrobora con la porosidad de 
retención de humedad, ya que todas las 
mezclas están entre 60% y 80%, aproxi- 
madamente (Tabla 3). Esto significa que 
hay muchos microporos llenos de agua 
que no estarían disponibles a las raíces 
de la planta. 

En el medio de crecimiento, el agua 
tiene un comportamiento que está en 
función de los sustratos usados para 
preparar las mezclas, así la disponibilidad 
del agua determinará si la mezcla será 
adecuada para un determinado cultivo. 
Las mezclas indican variaciones en los 
diferentes componentes que integran la 
curva de liberación de agua (Tabla 4). 

En la figura 1 se muestran las dife- 
rentes proporciones de los componentes 
de las curvas de liberación de agua para 
las cinco mezclas analizadas. El agua 
difícilmente disponible (add) en la mezcla 
de 100% de corteza es la más alta y 
disminuye al aumentar los porcentajes 
de aserrín y al bajar los porcentajes de 
corteza, lo que indica que mayores 
porcentajes de corteza en una mezcla 
tienden a retener el agua. Por tanto, el 
agua fácilmente disponible (afd) en las 
mezclas que contienen corteza es de las 
más bajas, con mayor materia sólida (ms) 
y que disminuye conforme se reducen los 
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Tabla 4. Componentes de la curva de liberación de agua en seis mezclas 
a base de corteza y aserrín, turba, agrolita y vermiculita. 


Mezcla Porcentaje de componentes 
de las mezclas 

ADD 

Componentes de la curva (%) 
AR AFD CA EPT 

MS 

1 

Corteza 100% 

38,8 

5,6 

6,8 

18,8 

70 

30 

2 

Corteza 80% + Aserrín 20% 

34,1 

5,1 

4,8 

33,0 

77 

23 

3 

Corteza 60% + Aserrín 40% 

35,9 

6,1 

7,3 

28,7 

78 

22 

4 

Corteza 20% + Aserrín 80% 

25,7 

4,2 

9,9 

51,2 

91 

9 

5* 

Turba 60% + Agrolita 30% + Vermiculita 10% 

26,0 

10,3 

17,7 

30,0 

84 

16 


*= Mezcla testigo, add= Agua difícilmente disponible, ar= Agua de reserva, afd= Agua fácilmente disponible, 
CA= Capacidad de aire, ept= Espacio poroso total y ms= material sólido. 


porcentajes de corteza. El afd en todas 
las mezclas está por debajo de 20%, que 
es lo mínimo que sugieren Abad et al. 
(1993) para mantener una buena propor- 
ción de agua en el sustrato. El agua de 
reserva (ar) está entre 4% y 10%, que es 
lo que sugieren Abad et al. (1993). 
Sumando el afd y el ar de la mezcla 
testigo, se llega a obtener 28% de agua 
total disponible. Esto es, un tercio del 
volumen del espacio poroso total (ept), 
más de un tercio es de capacidad de aire 
(ca) y un tercio aproximadamente es de 

ADD. 

CONCLUSIONES 

De acuerdo a las propiedades 
químicas, la mezcla dos (80% de corteza 
+ 20% de aserrín) es la que más se 
acerca a la mezcla estándar (60% de 
turba + 30% de agrolita + 10% de vermi- 
culita) en su pH, conductividad eléctrica y 
capacidad de intercambio catiónico. Sin 
embargo, la mezcla cuatro (20% de 
corteza + 80% de aserrín) presenta 
propiedades físicas como porosidad total, 
porosidad de aireación, porosidad de 
retención de humedad y agua fácilmente 


disponible, caacterísticas semejantes a la 
mezcla testigo. La mezcla cuatro 
presenta características prometedoras 
para ser utilizada en la producción de 
plantas forestales en contenedores 
rígidos (137 cm^) en vivero. 
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Estructura de bosques de pino pátula 
bajo manejo en Ixtlán de Juárez, 
Oaxaca, México 

Forest structure of managed patula pine in Ixtlán de Juárez, 

Oaxaca, México 
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Garza2, Óscar Alberto Aguirre-Calderón2, Javier Jiménez-Pérez2, 
Miguel Musalem-Santiago'i y Ricardo López-Aguillón2 

RESUMEN 


Se evaluó y comparó la diversidad de especies arbóreas y la composición estructural de 
bosques de pino pátula {Pinus patula SchI et Cham) bajo diferentes condiciones silvícolas, en Ixtlán 
de Juárez, Oaxaca, México. Se registró información sólo para ejemplares mayores de 7,5 cm de 
diámetro normal. La diversidad se estimó a través del Índice Shannon (H’), la composición estructural 
se definió mediante tres características, mezcla de especies, distribución espacial y diferenciación de 
tamaños, a través de los Indices, Mezcla de Gadow (Mi), Uniformidad de Ángulos (Wi) y Dominancia 
(Ui), respectivamente. El análisis se realizó comparando cuatro condiciones silvícolas denominadas, 
latizal, fustal joven, fustal medio y fustal viejo. Los resultados mostraron que la mayor diversidad 
arbórea se presentó en los bosques fustal medio y fustal viejo. El mayor grado de mezcla se presentó 
en la condición fustal viejo. La distribución del arbolado en el terreno, en los cuatro ambientes, 
presentó un patrón de distribución aleatorio y la dominancia de pino pátula va en aumento conforme 
la condición se desarrolla. 

PALABRAS CLAVE: 

Diversidad arbórea, índices estructurales, Pinus patula SchI ef Cham. 


ABSTRACT 


The diversity of arboreal species and the forest structure composition of Pinus patula SchI et 
Cham under different silvicultura! conditions, in Ixtlán of Juárez, Oaxaca, México was evaluated and 
compared. The information was registered oniy for trees over 7,5 cm of diameter breast height (dbh). 
The diversity was estimated through the Shannon Index (H’), the structural composition was defined 
through three characteristics, mingling species, spatial distribution and size differentiation, through the 
indices, Mingling by Godow (Mi), Contagien (Wi) and Dominance (Ui), respectively. The analysis was 
carried out comparing four silvicultura! conditions called, pole stage, young stem, médium stem and 
oíd stem. The comparative results showed that greater diversity was presentad in the médium stem 
and oíd stem forest. The major mingling degree was presentad in the oíd stem condition. The tree 
distribution, in the four environments, presentad a random distribution pattern and the dominance of 
Pinus patula is increasing while the condition develops. 

KEY WORDS: 

Pinus patula SchI et Cham, species diversity, structural indices. 
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INTRODUCCIÓN 

El estado de Oaxaca, México, es 
considerado uno de los centros de diver- 
sidad biológica más importantes a nivel 
nacional (García et al., 2004). En parti- 
cular para la región norte del estado, en 
la que se ubica Ixtlán de Juárez, los 
ecosistemas de pino-encino son de gran 
importancia ecológica, económica y 
social. En esta región, la mayoría de los 
bosques de pino-encino se conforman de 
rodales de segundo crecimiento que han 
estado bajo aprovechamiento comercial 
a lo largo de seis décadas aproximada- 
mente. Los principales tratamientos silví- 
colas han consistido en la aplicación de 
cortas selectivas y/o aclareos con el 
objetivo de producir madera aserrada y 
celulosa (uief Fapatux, 1977). A partir de 
1993 se incorporaron cortas de regenera- 
ción con el método de árboles padres 
(TIASA, 1993). 

En Ixtlán de Juárez los bosques de 
pino pátula {Pinus patula SchI et Cham) 
son los de mayor extensión, se distri- 
buyen en aproximadamente 5000 ha y es 
la especie de mayor importancia econó- 
mica. La medición cuantitativa de 
aspectos relacionados con la estructura 
de los rodales es muy importante para 
entender el funcionamiento del ecosis- 
tema, lo cual puede aportar elementos de 
decisión para contribuir al manejo 
forestal. 

La estructura de la masa forestal se 
puede definir como la forma en que los 
diferentes elementos del sistema se orga- 
nizan en el espacio (Kimmins, 1997) y se 
puede describir mediante tres caracterís- 
ticas: diversidad y mezcla de especies; 
distribución espacial y variaciones en las 
dimensiones de los árboles (Gadow y 
Huí, 1999; Frankiin et al., 2002). 

En la característica de diversidad 
interesan tres aspectos, el número de 


especies, la proporción de cada una de 
ellas, así como el grado de mezcla o la 
forma en que las especies se distri- 
buyen en el espacio. La distribución 
espacial está condicionada, principal- 
mente, por las relaciones entre indivi- 
duos, las estrategias de regeneración 
de las diferentes especies, las interven- 
ciones silvícolas y la fase de desarrollo 
del rodal. La variación en las dimen- 
siones de los árboles, describe la rela- 
ción de tamaños entre árboles vecinos a 
través de distintas variables, siendo las 
más frecuentes el diámetro y la altura 
(Del Río et al., 2003). 

En las áreas bajo manejo, los trata- 
mientos silvícolas modifican la estructura 
de los bosques, es por ello que su cono- 
cimiento, aunado a la dinámica de 
rodales, es fundamental para garantizar 
la sostenibilidad de estos sistemas. En 
este sentido, se planteó como hipótesis 
la existencia de diferencias significativas 
en la diversidad de especies arbóreas y 
en la composición estructural de los 
bosques de pino pátula, de acuerdo a la 
condición silvícola en que se presentan 
los rodales. 

Se ha definido que en bosques 
manejados, la biodiversidad es un 
elemento clave para evaluar la estabilidad 
del sistema (Láhde et al., 1999). En 
algunos estudios (Crow et al., 2002; Ishii 
et al., 2004) se ha concluido que la simpli- 
ficación de la estructura disminuye las 
funciones del ecosistema y la complejidad 
estructural aumenta la productividad y 
biodiversidad. Aunque no siempre estruc- 
turas complejas implican mayor diver- 
sidad, existiendo ejemplos contrarios 
(Hunter, 1999). En México, los estudios 
sobre metodologías (Jiménez et al., 2001 ; 
Aguirre et al., 2003) y evaluaciones del 
efecto que causan los tratamientos silví- 
colas (Corral et al., 2005; Solís et al., 
2006) sobre la diversidad estructural son 
todavía escasos. 
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OBJETIVO 

Evaluar y comparar la diversidad de 
especies arbóreas y la composición 
estructural de bosques de pino pátula 
{Pinus patula SchI et Cham), bajo dife- 
rentes condiciones silvícolas, en Ixtlán de 
Juárez, Oaxaca, México. 

METODOLOGÍA 

Ubicación 

El estudio se realizó en rodales de 
pino-encino con dominancia de pino 
pátula, en el predio de Ixtlán de Juárez, 
Oaxaca, México, a 170 18’ 16” y l/o 34’ 
00” latitud norte; 96° 31’ 38” y 960 20’ 00” 
longitud oeste. En la región, el pino pátula 
se distribuye principalmente en un rango 
altitudinal entre 2000 y 2500 metros, bajo 
un clima templado húmedo y subhúmedo 
con lluvias en verano, con una tempera- 
tura promedio anual de 200C y una preci- 
pitación promedio anual de 1100 mm. 

Diversidad de especies arbóreas 

En el bosque de pino pátula se imple- 
mentó un diseño de muestreo al azar 
estratificado, en donde los estratos fueron 
rodales en diferentes condiciones de 
desarrollo silvícola, definidos como: latizal, 
fustal joven, fustal medio y fustal viejo. 

Los criterios de selección de cada 
condición silvícola consistieron en la 
dominancia de pino pátula, que tuvieran 
una buena cobertura, con el mayor 
número de árboles ubicados en las clases 
diamétricas de 20 cm, 35 cm, 40 cm y 55 
cm, para cada condición, respectivamente, 
y ser bosques bajo aprovechamiento en 
diferentes tiempos, con la aplicación prin- 
cipalmente de cortas selectivas. 

El muestreo cubrió un área total de 
22,000 m2. En cada condición silvícola se 
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registró información en 5500 m2, resul- 
tado del levantamiento de 11 sitios de 
muestreo, cada sitio con un área de 500 
m2, de forma circular. Se registró informa- 
ción solo para las especies arbóreas 
mayores de 7,5 cm de diámetro, las varia- 
bles cuantificadas fueron: la especie, el 
número de árboles, el diámetro normal y 
la altura total, así como las distancias 
y ángulos del árbol central a cada uno de 
los árboles dentro del sitio. 

A partir de la información anterior, 
se derivó información fitosociológica y 
dasométrica, así para cada una de las 
especies se cuantificó su dominancia a 
través del área basal, su abundancia de 
acuerdo al número de árboles y su 
frecuencia a través de la presencia en las 
parcelas de muestreo; con estos tres 
datos se generó un valor ponderado para 
cada especie, denominado Valor de 
Importancia (vi). 

La diversidad de especies arbóreas 
se estimó a través del índice de Shannon 
(Magurran, 2004), el cual es un algoritmo 
que describe la diversidad de especies 
independiente de la distancia. El índice 
refleja la heterogeneidad de una comu- 
nidad sobre la base de dos factores: el 
número de especies presentes y su abun- 
dancia relativa. La fórmula correspon- 
diente aparece en la tabla 1 . 

Composición de la estructura 

Las estructuras de las condiciones 
silvícolas, latizal, fustal joven, fustal 
medio y fustal viejo, se describieron 
mediante tres características: mezcla de 
especies, distribución espacial y diferen- 
ciación de tamaños. El fundamento para 
el cálculo de los índices descriptivos, 
consistió en utilizar como unidad de 
muestreo al grupo estructural, que es un 
conjunto de cinco árboles, de los cuales 
un árbol sirve de referencia e incluye los 
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cuatro árboles vecinos más cercanos con 
los que convive (Pommerening, 2002; 
Kint et al., 2003). 

El tamaño de muestra para este 
análisis lo constituyeron 1176 grupos 
estructurales, los cuales fueron derivados 
de la información de distancias y ángulos 
contenidos en los 44 sitios de muestreo 
(11 por condición silvícola), a través de su 
ubicación por medio de coordenadas 
cartesianas y la aplicación de funciones 
trigonométricas. Se excluyeron del 
análisis aquellos árboles de referencia 
ubicados cerca del límite del sitio, cuyos 
árboles vecinos se encontraban fuera del 
radio de interés y con ello reducir la esti- 
mación del error (Pommerening y Stoyan, 
2006). Para describir el grado de mezcla 
de las especies arbóreas se utilizó el 
índice Mezcla de Gadow (1993), el cual 
se define como la probabilidad de que 
alguno de los cuatro árboles vecinos más 


cercanos sea de la misma especie que el 
árbol de referencia. 

Para cuantificar la distribución espa- 
cial horizontal de los árboles se utilizó el 
índice Uniformidad de Ángulos (Gadow y 
Huí, 2002), el cual se basa en la medición 
de ángulos entre dos vecinos (a) al árbol 
de referencia y considerando que la distri- 
bución horizontal de los árboles se ha 
clasificado típicamente bajo tres patrones 
generales: regular, aleatorio y en grupos 
(Zenner y Hibbs, 2000). 

Para cuantificar la variación exis- 
tente entre el tamaño de los árboles, en 
particular para la especie de pino pátula, 
y sus vecinos inmediatos, se utilizó el 
índice de Dominancia (Aguirre et al., 
2003). Para el cálculo se consideró el 
diámetro normal como variable cuantita- 
tiva. En la tabla 1 aparecen las fórmulas 
de los índices desarrollados. 


Tabla 1. Fórmulas de los índices de diversidad y estructurales utilizadas. 


índice 

Fórmula 

Significado 

Shannon (H’) 

H'=-Ípj*(lnpj) 

i = 1 

Pj = abundancia relativa de cada 
especie; n = número de especies. 

Mezcla de Gadow (Mi) 

1 N 14 

Mi = -2 -2 Vil 

Ni 4 J=1 

vy = 0 cuando el árbol vecino j perte- 
nece a la misma especie que /; 1 


de otra forma. 

Uniformidad de Ángulos (Wi) 

1 N 1 n 

Wi = -2 -2^1 

Np np 

Wj = 1 cuando el ángulo a es más 
pequeño que 720; o de otra forma. 

Dominancia (Ui) 

1 4 

Ui = -2vi 

4j=i 

Vj = dominancia del árbol vecino j, su 
valor es 1 cuando el árbol es más 



pequeño en diámetro que el árbol 
de referencia /; 0 de otra forma. 
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Análisis estadístico 

Con la finalidad de comparar esta- 
dísticamente los índices de diversidad y 
de estructura entre las condiciones silví- 
colas identificadas, se amplió la base de 
datos muéstrales, a partir de datos ya 
existentes, por medio de la técnica de 
remuestreo “bootstrap" estándar (Efron, 
1979). A partir de la información de 1000 
submuestras, se determinó la distribución 
de probabilidades de los índices a través 
de la prueba de Kolmogorov-Smirnov 
(ks). Con base en la prueba ks, se decidió 
llevar a cabo un análisis de varianza no 
paramétrico de un factor con la prueba de 
Kruskal Wallis, bajo la hipótesis nula de 
igualdad de medianas y se complementó 
con la prueba de comparaciones múlti- 
ples de Nemenyi (Zar, 1999). Para ambas 
pruebas se planteó un nivel de signifi- 
cancia de 95%. 

RESULTADOS Y DISCUSIÓN 
Diversidad de especies arbóreas 

Se registraron para toda la masa 
forestal 26 especies. De este total, en el 
bosque latizal se registraron 11 especies; 
en el bosque fustal joven se cuantificaron 
15 especies; fustal medio, 17 especies y 
en el bosque fustal viejo, 1 3 especies. De 
las 26 especies, sólo 7 de ellas estuvieron 
presentes de forma común, es decir, en 
las cuatro condiciones, y fueron las 
mismas que mostraron los mayores vi 
(Tabla 2). Destaca pino pátula como la 
especie más importante, lo cual significa 
mayor dominancia estructural y un nivel 
de competitividad mayor dentro de la 
comunidad. 

Las 19 especies restantes y con 
presencia alterna en las cuatro condi- 
ciones presentaron vi menores de 3%, 
esto significa que, a pesar de existir un 
número razonable de especies, en la 
mayoría de los casos se trata de árboles 
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en cantidades mínimas, de pequeñas 
dimensiones y distribuidos esporádica- 
mente. 

Para la condición latizal, el número 
total de árboles fue de 780 y el área 
basal (ab) de 25,3 m2ha-i; en la condi- 
ción fustal joven se registraron 607 
árboles y un ab de 37,6 m2ha-‘'; la 
condición fustal medio presentó 665 
árboles y el ab de 47,4 m2ha-’i y en la 
condición fustal viejo el número de 
árboles fue de 597 y el ab de 50 m2ha-L 
En la figura 1 se presenta la distribución 
diamétrica de pino pátula en cada una 
de las condiciones. En la tabla 2 se 
observa que el número de árboles para 
las especies de pino decrece conforme 
la condición fitosociológica se desarrolla, 
ocasionado por la mortandad de árboles, 
sobre todo de categorías pequeñas, atri- 
buible a la alta competencia y a la intole- 
rancia a la sombra. Para el caso de las 
latifoliadas, se mantiene e inclusive en 
algunas especies se presenta una 
tendencia a la alza, debido a la capa- 
cidad de rebrote que presentan, al incre- 
mento en diámetro y su incorporación en 
categorías mayores. 

Los valores de la diversidad de 
especies, cuantificada a través de H’ se 
presenta en la tabla 3. La prueba de 
Kolmogorov-Smirnov mostró que H’ no se 
distribuye normalmente (ks = 0,0438 
p = 0,0001). La prueba de Kruskal-Wailis 
mostró diferencias estadísticas significa- 
tivas en términos de diversidad arbórea 
entre las condiciones (chi-cuadrada 
= 251 8,39 y un valor p = 0,0001 ). A través 
de la prueba de Nemenyi se definió que la 
mayor diversidad arbórea se presentó en 
los bosques fustal medio y fustal viejo. 

A pesar que en la condición fustal 
medio se registró mayor número de espe- 
cies que en la condición fustal viejo, los 
valores de H’ estadísticamente perte- 
necen al mismo grupo de significancia. 
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Tabla 2. Valor de importancia (vi), número de árboles (na), área basal (ab), altura (a) 
y diámetro (di 3 m) para las especies comunes en cada condición silvícola 


ESPECIE 

VI 

NA.ha-1 

AB 

Altura (m) 

Diámetro (cm) 


% 


m2.ha*i 

A 

coef.var 

di.3 m 

coef.var 

Condición latizal (I): 








Pinus patula SchI et Cham. 

51 

490 

17,0 

16 

36 

21 

50 

Quercus crassifolia Bonpl. 

11 

73 

3,6 

14 

22 

25 

22 

Quercus rugosa Masam. 

8 

58 

1,5 

10 

25 

18 

32 

Pinus ayacahuite Ehrenberg. 

12 

100 

2,0 

16 

23 

16 

32 

Pinus pseudostrobus Lindiey. 

7 

20 

0,8 

17 

23 

22 

32 

Quercus corrúgate Hook. 

2 

10 

0,2 

8 

40 

17 

34 

Clethra tanate Michx. 

4 

15 

0,2 

9 

25 

13 

9 

Condición fustal joven (II): 








Pinus patula SchI eíCham. 

40 

236 

21,4 

23 

27 

34 

40 

Quercus crassifolia Bonpl. 

16 

87 

5,3 

12 

13 

28 

11 

Quercus rugosa Masam. 

11 

89 

3,1 

14 

33 

21 

25 

Pinus ayacahuite Ehrenberg. 

10 

73 

3,0 

16 

26 

23 

30 

Pinus pseudostrobus Lindiey. 

3 

18 

0,7 

15 

41 

22 

45 

Quercus corrúgate Hook. 

5 

27 

2,0 

15 

24 

31 

66 

Clethra ¡anata Michx. 

4 

27 

0,4 

9 

10 

13 

12 

Condición fustal medio (III): 








Pinus patula SchI et Cham. 

39 

233 

29 

26 

16 

40 

26 

Quercus crassifolia Bonpl. 

10 

98 

4 

11 

35 

23 

49 

Quercus rugosa Masam. 

6 

75 

1 

12 

31 

15 

46 

Pinus ayacahuite Ehrenberg. 

6 

63 

3 

17 

20 

24 

25 

Pinus pseudostrobus Lindiey. 

6 

18 

2 

28 

21 

39 

35 

Quercus corrugata Hook. 

8 

38 

1 

12 

25 

18 

22 

Clethra ¡anata Michx. 

2 

25 

0 

8 

20 

13 

10 

Condición fustal viejo (IV): 








Pinus patula SchI et Cham. 

32 

123 

28,2 

34 

15 

54 

25 

Quercus crassifolia Bonpl. 

15 

123 

7,1 

15 

26 

27 

31 

Quercus rugosa Masam. 

13 

110 

4,6 

16 

31 

23 

37 

Pinus ayacahuite Ehrenberg. 

7 

43 

2,1 

18 

20 

25 

35 

Pinus pseudostrobus Lindiey. 

7 

17 

4,4 

35 

8 

58 

19 

Quercus corrugata Hook. 

6 

23 

0,5 

14 

13 

16 

14 

Clethra ¡anata Michx. 

4 

23 

0,4 

11 

32 

14 

16 


otras especies y ubicación: Cinnamomum areolatum (L.) Kosterm (II, III, IV); Psychotria sp. (I, II, III); Arbutus xala 
pensis Kunth (I, II, III); AInus acuminata Kunth (III, IV); Pinus teocote SchI et Cham (III); Dendropanax arboreus (L.) 
Decne et lanch (IV); Pinus douglasiana Mtz (II, III); Cinnamomum sp. (III, IV), Pinus oocarpa Schiede (III); Abies 
hickeli Flous et Gaussen (IV); Quercus laurina Bonpl (II); Litsea glaucescens H.B.K.(I); Litsea sp. (III); Prunus sero 
tina Ehrh (II); Pinus oaxacana Mirov (I); Ternstroemia pringlei Rose (IV); Quercus scytophylla Liebm (II); Rheedia 
macrophylla (Mtz.) Planch ef Triana (II); Viburnum sp, (III). 
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Categoría diamétrica (cm) 


— ♦ — Latizal 
— ís — Fustal joven 
- - - ■ Fustal medio 
— e — Fustal viejo 


Figura 1. Distribución de categorías diamétricas de pino pátula por condición. 


atribuidles a las mayores abundancias 
relativas de las especies latifoliadas de la 
última condición. Una situación similar se 
presentó entre las condiciones fustal 
joven y fustal viejo. Lo anterior confirma 
que H’ aumenta de valor conforme es 
mayor el número de especies presentes y 
también cuando las proporciones de las 
distintas especies son similares. 

De lo anterior se deriva que la apli- 
cación de tratamientos silvícolas ha modi- 
ficado las condiciones de diversidad y 


abundancia de las especies arbóreas 
presentes. En bosques templados de 
Durango, Solís et al. (2006) registraron 
también que las cortas selectivas dismi- 
nuyeron la diversidad arbórea. 

Composición de la estructura 

En lo que corresponde al grado 
de mezcla de las especies, la prueba de 
Kruskal-Wailis detectó la existencia de 
diferencias estadísticas significativas (chi- 


Tabla 3. Valores promedio del índice Shannon (H’), rango medio de la prueba 
de Kruskal-Wailis (kw) y comparación múltiple con la prueba de Nemenyi 

por condición silvícola. 


Condición 

silvícola 

Número de 
Especies 

índice 

H’ 

Rango Medio 
de KW para H’ 

Prueba de 
Nemenyi 

Latizal 

11 

1,3 

631 

c 

Fustal joven 

15 

1,9 

1716 

b 

Fustal medio 

17 

2,17 

2965 

a 

Fustal vieio 

13 

2,16 

2687 

a 
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Tabla 4. Valores promedio para los índices, Mezcla de Gadow (Mi), Uniformidad 
de Ángulos (Wi), Dominancia (Ui) y rango medio de la prueba de Kruskal-Wailis (kw) 
y comparación múltiple con la prueba de Nemenyi por condición silvícola. 


Condición 

índice Mi 

Rango Medio 
de KW y Prueba de 
Nemenyi para Mi 

índice Wi 

Rango Medio 
de KW y Prueba de 
Nemenyi para Wi 

índice Ui 

Rango Medio 
de KW y Prueba de 
de Nemenyi para Ui 

Latizal 

0,45 

344c 

0,50 

401a 

0,52 

159c 

Fustal joven 

0,56 

425b 

0,51 

414a 

0,56 

171c 

Fustal medio 

0,58 

432b 

0,54 

434a 

0,70 

208b 

Fustal vieio 

0,69 

514a 

0,52 

419a 

0,83 

249a 


cuadrado = 53.7412 y p = 0.0001) entre 
las condiciones bajo estudio (Tabla 4). A 
través de la prueba de Nemenyi se definió 
que, la condición fustal viejo con un índice 
Mi de 0,70, se mezcla en mayor grado 
que el resto de las condiciones. Ello signi- 
fica que en promedio, de los cuatro 
árboles circunvecinos más cercanos al 
árbol de referencia, tres de ellos perte- 
necen a especies diferentes. 

El índice de mezcla está determi- 
nado principalmente por la abundancia 
relativa de las especies de árboles. Espe- 
cies presentes en alta proporción refle- 


jarán valores bajos de mezcla, como es 
de suponer al tener pocos vecinos de 
especies diferentes. 

Para un análisis con mayor detalle 
de la mezcla espacial de los árboles, se 
representó gráficamente la distribución de 
los valores del índice Mi, para cada una 
de las condiciones silvícolas (Figura 2). 
En la condición latizal se observa una 
tendencia a mantener una población 
menos mezclada. El 23% del arbolado 
tiene un valor de mezcla igual a cero, lo 
que significa que estos árboles 
conforman grupos de cinco individuos de 



□ Latizal 
e Fustal joven 
ii Fustal medio 
■ Fustal viejo 


Figura 2. Distribución de la mezcla de especies por condición silvícola. 
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una misma especie y 35% del arbolado 
se mezcla con tres y cuatro especies dife- 
rentes. En contraste, la condición fustal 
viejo presenta 2% y 61%, respectiva- 
mente. 

En lo que respecta al patrón de 
distribución horizontal de los árboles, los 
valores promedio del índice Wi a nivel de 
condición se presentan en la tabla 4. La 
prueba de Kruskal-Wailis indicó que no 
existen diferencias estadísticas entre 
condiciones (chi-cuadrado = 2,8608 y 
p = 0,4136). Los cuatro ambientes pre’- 
sentaron un patrón de distribución alea- 
torio. De acuerdo con Hui y Gadow 
(2002), el valor medio del índice Wi de 
alrededor de 0,50, significa que el arreglo 
espacial horizontal de los árboles pre’- 
senta un patrón general de distribución 
aleatorio, conforme el valor aumenta de 
0,6 la distribución tiende a ser en grupos 
y en un valor menor de 0,4 la distribución 
es regular, como en plantaciones jóvenes. 

Los valores del índice Wi aquí repor- 
tados coinciden con otros estudios como 
el de Aguirre et al. (2003), quienes en tres 
bosques de coniferas excluidos del 
manejo forestal, definieron la distribución 
de los árboles como aleatoria, al deter- 


minar valores promedio para el índice Wi 
de 0,50, 0,52 y 0,51. Pommerening 
(2002) reportó en tres plantaciones 
adultas un patrón de distribución alea- 
torio, con valores Wi de 0,57, 0,53 y 0,56, 
respectivamente, atribuidle a la 
mortandad, establecimiento natural de 
otras especies y actividades de manejo 
forestal. 

La distribución del índice Wi para las 
cuatro condiciones silvícolas puede obser- 
varse en la figura 3. En ella se muestra una 
tendencia similar con pequeñas varia- 
ciones. Para las cuatro condiciones, en 
promedio, el porcentaje de árboles con 
una distribución regular fue de 21%, con 
distribución aleatoria 54% y con una distri- 
bución en grupos de 25%. Al analizar la 
distribución espacial horizontal de los 
árboles solamente de pino pátula, ésta se 
mantuvo como aleatoria en las cuatro 
condiciones, estadísticamente tampoco se 
detectaron diferencias significativas (Chi- 
cuadrada = 3.2440 y p > 0.3555). 

Investigaciones afines (Ishii et al., 
2004; Corral et al., 2005) han registrado 
cambios en la distribución espacial hori- 
zontal de los árboles, de una distribución 
regular en áreas sin aprovechamiento a 



□ Latizal 
ri Fustal joven 
u Fustal medio 
■ Fustal viejo 


Figura 3. Distribución del índice de Uniformidad de Ángulos (Wi) por condición silvícola. 
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□ Latizal 
s Fustal joven 
tt Fustal medio 
■ Fustal viejo 


Figura 4. Distribución del índice de Dominancia de pino pátula por condición silvícola. 


una distribución en grupos en áreas bajo 
aprovechamiento forestal. Al respecto, 
Neumann (2001) menciona que la distri- 
bución regular proporciona el máximo 
espacio y condiciones de crecimiento 
óptimas para cada individuo, mientras la 
distribución en grupos provoca pérdidas 
de incremento. Los patrones de distribu- 
ción han sido vinculados a diferentes 
procesos (como la mortalidad de árboles, 
la competencia, regeneraciónj creación de 
huecos, seminación, entre otros) de ser 
los responsables de los patrones espa- 
ciales observados (Zenner y Hibbs, 2000). 

En lo que corresponde a la diferen- 
ciación de tamaños en diámetro para la 
especie de pino pátula cuantificada por 
medio del índice Di, la prueba de Kruskal- 
Wallis detectó la existencia de diferencias 
estadísticas significativas (chi-cuadrado = 
25.1586 y p = 0.0001) entre las condi- 
ciones silvícolas bajo estudio. En la tabla 
4 se registran los valores promedio de Di 
para cada condición, se observa que 
conforme la condición se desarrolla, la 
dominancia de pino pátula va en 


aumento. Para las condiciones latizal y 
fustal joven, pino pátula es codominante, 
en promedio, se asocia con dos árboles 
más delgados y con dos más gruesos. En 
la condición fustal medio, pino pátula es 
dominante al presentarse entre dos y tres 
árboles circunvecinos más delgados. En 
la condición fustal viejo, pino pátula es 
ampliamente dominante ya que al menos 
tres árboles de los cuatro que lo rodean 
son más delgados. 

En la figura 4 se representa con 
mayor detalle la distribución del índice Di 
en cada una de las cuatro condiciones. 
En la condición latizal, la distribución 
porcentual de pino pátula se mantuvo de 
una forma equilibrada, con pequeñas 
variaciones alrededor de 20% en cada 
uno de los cinco criterios de clasificación. 
En la condición fustal viejo, pino pátula es 
una especie plenamente dominante, ya 
que 29% de los árboles de referencia 
están rodeados de tres árboles más 
pequeños y 54% de cuatro árboles más 
pequeños. 
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CONCLUSIONES 

La cuantificación de aspectos rela- 
cionados con la estructura de los rodales 
es importante para entender el funciona- 
miento del ecosistema y con ello contri- 
buir al manejo forestal sustentadle. Se 
presentaron diferencias en la diversidad 
de especies arbóreas según la etapa de 
desarrollo silvícola. En bosques maduros, 
del tipo fustal medio y fustal viejo, se 
registraron valores que indicaron mayor 
diversidad de árboles que en bosques 
inmaduros del tipo latizal y fustal joven, 
que mostraron menor diversidad. Pino 
pátula es la especie más importante en 
cada una de las condiciones silvícolas, al 
ocupar mayor área basal, registrar el 
mayor número de árboles por unidad de 
superficie y de mayor distribución en el 
área de estudio. El grado de mezcla de 
las especies arbóreas varió de acuerdo a 
la condición silvícola, éste avanzó progre- 
sivamente conforme la condición se 
desarrolló. En el tipo de bosque latizal se 
presentó la condición más pura o el 
menor grado de mezcla, hasta el bosque 
de tipo fustal viejo donde se presentó la 
mayor mezcla de las especies. 

El arreglo espacial horizontal de los 
árboles se comportó de manera similar en 
las cuatro condiciones silvícolas. Inde- 
pendientemente del desarrollo de la 
masa, mostró un patrón general de distri- 
bución aleatorio, aunque también se 
registraron proporciones menores con 
distribución regular y en grupos. La domi- 
nancia en diámetro (y en consecuencia 
en área basal) de pino pátula sobre sus 
vecinos más cercanos varía y crece con- 
forme la condición silvícola se desarrolla, 
desde la condición latizal donde es codo- 
minante hasta la condición fustal viejo 
donde es plenamente dominante. Final- 
mente, cabe mencionar que la aplicación 
de tratamientos silvícolas ha modificado 
las condiciones de diversidad y estructura 
de rodales. 
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lARTICULO DE INVESTIGACION 


Crecimiento y rendimiento de Cedrela odorata L 
y Tabebuia donnell-smithii Rose 
en San José Chacaiapa, Pochutia, Oaxaca 

Growth and wood yieldíng of Cedrela odorata I. 
and Tabebuia donnell-smithii Rose ín San José Chacaiapa, 

Pochutia, Oaxaca 

Rolando Galán Larrea^, Héctor Manuel de los Santos Posadas2, 

y Juan Ignacio Valdez Hernández2 

RESUMEN 

Se construyó un Sistema de Crecimiento Maderable (SCM) de tipo explícito a partir de datos de 
sitios con remediciones de Cedrela odorata (cedro rojo) y Tabebuia donnell-smithii (primavera) en la 
región de la costa oaxaqueña. Como primera fase del SCM se construyeron modelos de volumen total 
y comercial por árbol con base en datos de volumen obtenidos mediante métodos no destructivos. Los 
datos de las parcelas se ajustaron a un SCM explícito, a través de regresión aparentemente no corre- 
lacionada para lograr compatibilidad total entre componentes de predicción y proyección. 

PALABRAS CLAVE; 

Diferencia algebraica, índice de Sitio, Pochutia. 


ABSTRACT 

An explicit growth and yieid system (GYS) was fitted to data from re-measured plots of Cedrela 
odorata (spanish cedar) and Tabebuia donnell-smithii plantad at the Oaxaca’s Coastal región. As a first 
stage, a total and merchantable volume system at tree level was fitted based on non-destructive 
sampling scheme. Afterwards a GYS was fitted using seemingly unrelated regression to achieve full 
compatibility among prediction and projection models. 

KEY WORDS: 

Algebraic difference. Site Index, Pochutia. 


INTRODUCCIÓN 

A pesar de existir una serie de inves- 
tigaciones encaminadas a generar 
nuevas técnicas para lograr un manejo 
óptimo de los recursos forestales, como 
los modelos de crecimiento y rendimiento, 


en México aún se tienen muchas caren- 
cias para el desarrollo, implementación y 
ejecución de estas herramientas, las 
cuales se han estudiado principalmente 
en especies arbóreas de bosques natu- 
rales en clima templado, como las coni- 
feras. 


1 Estudiante de Maestría, Posgrado Forestal (PF) Colegio de Posgraduados (COLPOS). ce: rgalan@colpos.mx 

2 PF-COLPOS. CP 56230, Montecillo, Texcoco, Estado de México. Tel. (595) 9520200 ext 1484, Tel. D.F. (55) 
58045900 ext 1484. ce: hmsantos@colpos.mx, ignaciov@colpos.mx 
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Para muchas especies arbóreas en 
los trópicos, no existe información cientí- 
fica previa acerca de su crecimiento y 
rendimiento, recurriendo comúnmente a 
estimaciones realizadas por los dueños del 
bosque. Esto lleva, a su vez, a una falta de 
información objetiva sobre el análisis 
benefido/costo para la relación especie- 
sitio de plantación y, como consecuencia, 
a la inexistencia de un verdadero manejo 
en la mayoría de los bosques tropicales 
(Ricker y Del Río, 2004). 

Las especies nativas de clima cálido 
en México, a pesar de su importancia 
económica y ecológica, carecen en su 
mayoría de información cuantitativa 
acerca de su crecimiento y desarrollo, 
básica para su manejo. Una alternativa 
para obtener dicha información es 
mediante el estudio de plantaciones esta- 
blecidas, en las cuales se conoce la edad 
de las especies plantadas. En México, se 
han realizado algunos trabajos relativos a 
sistemas de crecimiento y rendimiento 
maderable en plantaciones de Swietenia 
macrophylla (García et al., 1998), y 
Eucalyptus sp. (Gómez et al., 2006; 
Reyes et al., 2006). 

En la región de la costa de Oaxaca, 
fueron establecidas plantaciones fores- 
tales de enriquecimiento de acahual, con 
especies de alto valor comercial. Estas 
plantaciones siguen una distribución 
sistemática, regular y ordenada, en 
líneas, fajas o grupos dentro del bosque o 
acahual, y se consideran un método 
adecuado para regenerar bosques sobre- 
explotados donde la regeneración natural 
es insuficiente y se puede incrementar el 
volumen maderable con especies de alto 
valor comercial (Adjers et al., 1995; 
Montagníni et al., 1997). Sin embargo, 
actualmente no se cuenta con informa- 
ción sobre el crecimiento y rendimiento 
maderable de estas especies tropicales, 
sea en bosques naturales o en planta- 
ciones. 


OBJETIVO 

El objetivo del estudio es el 
desarrollo de un sistema de medición del 
crecimiento y rendimiento maderable 
para las especies Cedrela odorata (cedro) 
y Tabebuia donnell-smithii (primavera) en 
la región, a partir de datos de diámetro y 
altura obtenidos de remediciones en sitios 
de muestreo permanentes en planta- 
ciones de enriquecimiento de acahual. En 
una primera fase se determina un sistema 
de cubicación de volumen total del fuste 
para Cedrela odorata L. y Tabebuia 
donnell-smithii Rose, mediante modelos 
matemáticos. Posteriormente, se cons- 
truye un sistema de crecimiento y rendi- 
miento maderable, a través de modelos 
de predicción explicita (totalidad de rodal) 
para volumen por hectárea, y se evalúa la 
productividad de las especies estudiadas 
mediante la elaboración de un índice de 
sitio basado en el crecimiento de altura 
dominante. 

METODOLOGÍA 

La presente investigación se llevó a 
cabo en plantaciones de enriquecimiento 
de acahual con edades comprendidas de 
tres a 10 años, en los predios: El Pén’- 
jamo. El Carnizuelo, Arroyo Rico y El 
Riego, ubicados en las inmediaciones del 
poblado San José Chacalapa, pertene- 
ciente al municipio de Pochutia, estado 
de Oaxaca. 

El clima en el área de estudio, de 
acuerdo con el sistema de clasificación 
climática de Kóppen modificado por García 
1988, corresponde al grupo de climas 
cálidos subhúmedos con lluvias en verano 
Aw^: el intermedio entre Awq y Aw 2 , con 
una precipitación media anual de 934,7 
mm y una temperatura media anual de 
26°C (García, 1988). El suelo predomi- 
nante en el área de estudio corresponde a 
un regosol éutrico (inegi, 1995). 
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Actualmente el tipo de vegetación, 
en tres de los cuatro predios (El Pénjamo, 
Carnizuelo y Arroyo Rico), corresponde a 
acahuales formados como consecuencia 
de la roza, tumba y quema. Sin embargo, 
con base en la clasificación de la vegeta- 
ción de México realizada por Miranda y 
Hernández X. (1963), y tomando en 
cuenta las especies relicto, el tipo de 
vegetación fue una selva mediana caduci- 
folia, aproximadamente tres cuartas 
partes de los árboles altos pierden 
completamente sus hojas en la época de 
sequía (Miranda y Hernández X., op. cit). 

Para el caso del predio El Riego, 
todavía se observa la estructura tanto 
vertical como horizontal de la selva 
mediana caducifolia, aunque de mode- 
rada a alta la perturbación; el dosel de las 
especies vegetales presentes es utilizado 
como sombra de cafetales. 


Especies estudiadas 

Son nativas de México e importantes 
por sus valores ecológicos y comerciales; 
son descritas a continuación con base en 
el trabajo realizado por Pennington y 
Sarukhán (1998). 

Cedrela odorata L (MELIACEAE) 

Comúnmente se le conoce como 
cedro rojo, nombre aplicado en toda su 
área de distribución; aunque también es 
llamado como chujté (tzeltal) en la zona 
lacandona, Chiapas, así como icte (huas- 
teco) y tiocuáhuitl (náhuatl). Este árbol 
puede llegar a medir hasta 35 m de altura 
y alcanzar un diámetro normal de 1 ,7 m, 
en ocasiones forma pequeños contra- 
fuertes poco prominentes, sus ramas son 
ascendentes y gruesas, su copa es 
redonda y densa. Las hojas son 
dispuestas en espiral, paripinnadas o 
imparipinnadas, de 15 a 50 cm de largo 
incluyendo el pecíolo, compuestas por 10 
a 22 folíolos opuestos o alternos, además 


poseen un penetrante olor a ajo cuando 
se estrujan. Los árboles de esta especie 
tiran las hojas cuando han madurado 
totalmente los frutos de la temporada 
anterior, antes de florecer. 

Se encuentra en la vertiente del 
Golfo desde el sur de Tamaulipas y 
sureste de San Luis Potosí, hasta la 
península de Yucatán y en la del Pacífico 
desde Sinaloa hasta Guerrero, asi como 
en la Depresión Central y costa de 
Chiapas. Alcanza sus máximos tamaños 
e incrementos en zonas con precipita- 
ciones entre 2,500 y 4,000 mm anuales; 
en zonas con precipitación pluvial mucho 
menor, el árbol no se desarrolla tan bien y 
presenta fustes cortos y frecuentemente 
torcidos. 

El cedro después de la caoba es la 
especie maderable más importante en la 
Industria forestal de México. Su madera 
de características excelentes, se usa para 
obtener vigas, tablas y chapas, así como 
para fabricar diversos artículos tor- 
neados, cajas de puros y hacer tallas. El 
uso más indicado para esta especie, por 
el bello jaspeado de su madera, es la 
fabricación de chapas y madera terciada 
para exportación. 

Tabebuia donnell-smithii Rose (BIGNO - 
NIACEAE) 

Comúnmente se le conoce como 
primavera, en toda el área de distribución. 
Este árbol alcanza hasta 30 m de altura y 
un diámetro normal de 70 cm, con el 
tronco ligeramente acanalado; ramas 
ascendentes, copa alargada. Las hojas 
son decusadas digitado compuestas, de 
20 a 70 cm de largo incluyendo pecíolo, 
compuestas de siete folíolos, ocasional- 
mente cinco. Las flores se presentan en 
panículas terminales, piramidales, con 
ramas cimosas de 15 a 35 cm de largo. El 
fruto es una cápsula de 25 a 50 cm de 
largo. Los árboles de esta especie 
pierden las hojas a fines de marzo, época 



Madera y Bosques 14(2), 2008:65-82 

en que florecen, y las reponen a princi- 
pios de mayo. 

La especie se distribuye en la 
vertiente del Pacífico, desde Nayarit hasta 
Chiapas, y en la del Golfo en Tabasco, 
formando parte de las selvas medianas 
subcaducifolias en suelos de origen 
volcánico y metamórfico o sedimentario. 

Su madera, de color claro, es muy 
apreciada para la fabricación de 
muebles. 

Sitios permanentes de dimensiones 
fijas 

En marzo del 2005 se establecieron 
14 sitios permanentes de dimensiones 
fijas: 20 x 20 metros (400 m2), los cuales 
se ubicaron en el terreno utilizando cintas 
métricas, brújula (Suunto), machetes para 
limpia de límites y estacas de madera, las 
cuales fueron marcadas con pintura viní- 
lica. Los árboles ubicados dentro de los 
sitios fueron medidos en marzo del 2005 
y remedidos al año siguiente (marzo del 
2006). Cada árbol se identificó con 
pintura mediante un número secuenciado 
de control y se le registró su diámetro a la 
altura del pecho (dap), altura total y edad 
con base en los registros existentes de la 
plantación, asegurándose que correspon- 
dieran al paraje de estudio; además, se 
obtuvieron datos como nombre del 
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predio, número de sitio, especie, grado o 
porcentaje de inclinación de la pendiente, 
exposición de la pendiente y ubicación 
geográfica. 

Cubicación de árboles en pie 

Cedrela odorata se encuentra plan- 
tada en tres de los cuatro predios 
comprendidos para este estudio: El 
Pénjamo, El Carnizuelo y Arroyo Rico, en 
ellos se encontraron cinco, tres y dos 
categorías diamétricas de 5 cm, respecti- 
vamente; mientras que Tabebuia donnell- 
smithii solamente se encuentra estable- 
cida en el predio El Riego, con dos cate- 
gorías diamétricas de 5 cm. Se seleccio- 
naron ocho árboles por categoría diamé- 
trica y en cada predio (Tabla 1). 

Ante la falta de información referente 
a modelos de volumen de fuste para la 
zona, se midieron los diámetros normales 
con corteza, utilizando cinta diamétrica 
con un nivel de precisión de 3,14 mm 
(Forestry Suppliers, 1992), así como 
diámetros a diferentes alturas del fuste, 
iniciando en 0,30 cm y posteriormente 
cada metro. Los diámetros superiores 
fueron tomados de manera directa con 
ayuda de una escalera de madera de tres 
metros de alto, mientras que las alturas 
de fuste limpio y total del árbol se regis- 
traron con una baliza de madera 
graduada cada 50 cm. 


Tabla 1. Número de árboles de Cedrela odorata y Tabebuia donnell-smithii por categoría 
diamétrica, medidos en pie para obtener el sistema de cubicación. 


Especie 


Categoría diamétrica 

i ii iii iv V Totai 

(2,5-7,5cm) (7,6-12,5 cm) (12,6-17,5 cm) (17,6-22,5 cm) (22,6-27,5 cm) Totai 


Cedrela odorata 24 24 16 8 8 80 

Tabebuia donnell-smithii 8 8 16 
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La cubicación de árboles en pie es 
un procedimiento no destructivo, es decir, 
no es necesario derribar los árboles para 
cubicarlos, ya que se recurre a la medi- 
ción de alturas parciales y de diámetros 
en aquellos puntos seleccionados como 
divisiones entre troza y troza (Romahn et 
al., 1994). Se escogió este método por 
tratarse de especies con alto valor comer- 
cial y representar una de las primeras 
plantaciones establecidas en la zona. El 
volumen del fuste se estimó a partir de la 
cubicación de cada una de las trozas o 
segmentos en los que se dividió el árbol, 
sumando al final todos los segmentos; 
para la cubicación de cada segmento 
(troza) y tanto para volumen total como 
volumen comercial, se utilizó la fórmula 
de Smalian (Romahn, et al., op. cit.). 


Volumen total y comercial variable por 
árbol 


Para obtener el volumen comercial 
se ajustaron tres modelos matemáticos: 
Clutter, 1980 (2); Fang y Bailey, 1999 (3) 
y Piennar y Reheney, 1988 (4), los cuales 
involucran la estimación simultánea del 
volumen total y comercial variable. Estos 
modelos de volumen comercial ajustan 
simultáneamente con el modelo de 
volumen total de Schumacher y Hall que 
implica un modelo de ahusamiento. 


V=Kx 


l 






d^\ 


( 2 ) 


V^=V,-P„xd^x{H-h) (3) 


v. = v,-P,^ 


I d 


Pl \ 


D 


A-2 


X 


(íí-1.3) 


(4) 


Para obtener el volumen total con 
corteza, y tomando en cuenta la estruc- 
tura de los datos, se empleó el modelo no 
lineal de Schumacher y Hall (1933), ajus- 
tando el modelo de volumen total a través 
de las variables diámetro y altura, 
mediante análisis de regresión usando el 
paquete estadístico sas (Statistical 
Analysis System, 2001). Este modelo es 
una formulación que permite expresar el 
volumen de varios sólidos de revolución 
como el cono, paraboloide, neiloide y 
cilindro (Fang y Bailey, 1999); 


F, =aoxD“' (1) 


Donde Ve es el volumen comercial 
variable en m^, d es el diámetro en cm de 
punta comercial a la altura que h que es 
la altura comercial en m a partir del tocón; 
[5o, y P 2 corresponden a los parámetros 
a estimar para la punta no comercial, las 
demás variables fueron definidas en el 
modelo 1 . 

La elección del modelo que presentó 
mejor ajuste fue basada en los criterios 
tomados de Rentería (1995): alto valor de 
R2 o pseudo R2 (cercano a 1 ), valor bajo 
del Cuadrado Medio del Error, parsimonia 
y significancia estadística de los paráme- 
tros estimados. 


Donde: es el volumen total de fuste en 

m3, D es el diámetro normal o a la altura 
del pecho (dap) en cm, H es la altura en 
m, y son los parámetros a estimar. 


Estimación de alturas 

Para obtener la cubicación por sitio 
de muestreo es necesario medir el 
diámetro y altura total de todos los 
árboles presentes. Por razones de costo 
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y tiempo, así como en ocasiones por 
condiciones del terreno: alta inclinación 
de la pendiente, cárcavas y vegetación 
competidora (p. ej. malezas), esta medi- 
ción frecuentemente no se realiza. 

Ante tal situación, la altura total de 
una submuestra es utilizada en modelos 
de regresión para estimar las alturas del 
resto de la población en función del 
diámetro normal. La relación existente 
entre estas dos variables es muy 
estrecha, generalmente las alturas 
aumentan con el diámetro, sin embargo, 
llega un momento en que la altura se 
mantiene constante y el diámetro 
continua aumentando (García, 1995). En 
el presente estudio, durante la primera 
medición (marzo del 2005) fueron 
tomadas solamente las alturas de los 
árboles más altos (dominantes), esti- 
mando las alturas de aquellos árboles no 
medidos a través del siguiente modelo: 


D*’' (5) 

Donde: (po y (fi son los parámetros a 
estimar. 


Sistema de rendimiento maderable 

Los modelos de predicción explícita 
han sido ampliamente usados por su 
sencillez y precisión, requiriendo informa- 
ción básica como diámetro (dap) y altura 
para su construcción. Tienen la desven- 
taja de proporcionar información muy 
general de la plantación de interés; sin 
embargo, algunas modificaciones y 
ampliaciones de estos modelos pueden 
usarse para proporcionar información 
más detallada de la misma (Torres y 
Magaña, 2001). Por lo anterior, y conside- 
rando que en la región no se cuenta con 
información básica de estas especies, se 
utilizaron modelos de predicción explícita 
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para construir un sistema de crecimiento 
y rendimiento maderable (scrm), el cual 
permitirá generar datos acerca de la diná- 
mica del crecimiento en volumen de la 
plantación y como base para trabajos 
posteriores. 

Crecimiento en altura dominante 

A la altura dominante evaluada a 
una edad de referencia se le denomina 
índice de Sitio (is), y es el método directo 
más popular y práctico para la evaluación 
de la productividad forestal (Clutter et al., 
1983; Zepeda y Rivero, 1984). 

Para la construcción de las curvas 
de índice de sitio se utilizaron datos de 
altura y edad de árboles dominantes y 
codominantes en 14 sitios permanentes 
de muestreo, ubicados de tal modo que 
cubrieran los intervalos de edad (3-10 
años) y calidad de sitio. Primeramente se 
ajustó el modelo de la curva guía y poste- 
riormente se seleccionaron dos modelos 
a través de la estrategia denominada 
Método de la Diferencia Algebraica, final- 
mente se ajustaron simultáneamente los 
modelos de predicción y proyección, 
generando una familia de curvas de is de 
tipo anamórficas: la altura guarda la 
misma proporción a diferentes edades, 
haciendo que las curvas aparenten tener 
la misma forma (Clutter et al., 1983). 

Para la elaboración de estas curvas 
de índice de sitio, mediante el método de 
mínimos cuadrados, se ajustaron los 
modelos de Schumacher y Chapman- 
Richards, por ser sencillos y muy utili- 
zados al modelar altura dominante y 
construir este tipo de curvas (Ramírez, 
1981); además, mostraron buenos resul- 
tados de ajuste para modelar la altura 
dominante en Swietenia macrophylla 
mediante análisis de regresión (García et 
al., 1998). La selección del modelo más 
adecuado para describir la relación altura 
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dominante-edad se hizo mediante los crite- 
rios de pseudo R2, cuadrado medio del error 
(cme), valor de F, intervalos de confianza y 
error estándar de los estimadores. 

Las expresiones matemáticas, para 
los modelos promedio (6) y diferencia 
algebraica ( 7 ) muestran una familia de 
curvas anamórficas tipo Schumacher 
usada para estimar el is y la proyección 
de alturas en las plantaciones: 



Las expresiones para los modelos 
promedio (8) y diferencia algebraica ( 9 ) 
muestran una familia de curvas anamór- 
ficas tipo Chapman-Richards empleada 
para estimar el is y la proyección de 
alturas en las plantaciones: 

4 = Óq X (l - exp(- d^xE )J^ (8) 


del crecimiento de las especies en los 
índices de sitio definidos. 

Crecimiento en área basal 

El área basal es una medida directa 
de la densidad que involucra tanto el 
número de individuos como el tamaño de 
éstos, y es, quizá, la variable más utili- 
zada para modelar el efecto de la 
densidad en el crecimiento del arbolado, 
tanto a nivel poblacional como individual. 

Para obtener el área basal actual se 
utilizó el siguiente modelo no lineal 
(Torres y Magaña, 2001 ): 


-4 '|xA,* (10 ) 

E, j 


La proyección del área basal se 
realizó con el siguiente modelo: 



ABj = AB. X — 
A 


A + Él] ('11) 


Al — Al X 


1 - exp(- ó^xEi) 
1 - exp(- óixEi) 


Ó2 

( 9 ) 


Donde: A-/ es la altura dominante inicial 
en metros; A2 es la altura dominante a 
proyectar en metros; E es la edad en 
años; E7 es la edad inicial, E2 es la edad 
de proyección; óq, ói y 62 son los paráme- 
tros a estimar. 

En el presente estudio se seleccio- 
naron curvas de tipo anamórficas por 
presentar los mejores ajustes estadísticos 
y las mejores representaciones gráficas 


Donde ABi y A^ corresponden al área 
basal y altura dominante, respectiva- 
mente, a la edad inicial E7 (edad de refe- 
rencia); AB2 y A2 son área basal y altura 
dominante, respectivamente, a la edad E2 
(edad de proyección); y y 02 son los 
parámetros a estimar del modelo. 

Modelo de volumen total 

Para predecir y obtener el volumen 
total actual se empleó el siguiente modelo 
que utiliza el área basal y la edad: 

K=r<,x X expí “ ] (12) 
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También fue posible proyectar el 
volumen a través del modelo siguiente: 






compatibilidad total entre el modelo de 
estimación promedio y el de proyección, 
de forma que los parámetros estimados 
son los mismos mientras se cumple con 
el criterio de minimización de cuadrados. 


(13) 

Donde Vi y ABi son el volumen y área 
basal, respectivamente, a la edad inicial E-j 
(edad de referencia); \/2 y >AS 2 corresponden 
al volumen y área basal, respectivamente, a 
la edad E 2 (edad de proyección); y yo. Y 1 y Y 2 
son los parámetros a estimar. 


Método de ajuste del SCRM 

Antes de ajustar el sistema de creci- 
miento y rendimiento maderable, se partió 
del ajuste de modelos promedio para 
obtener buenos estimadores iniciales de 
los parámetros; posteriormente, se ajus- 
taron los modelos promedio y de proyec- 
ción como un sistema de ecuaciones bajo 
la técnica de regresión aparentemente no 
correlacionada (sur por sus siglas en 
inglés) (Borders, 1989). 

La técnica de estimación de paráme- 
tros SUR proporciona una mejor ganancia 
en la estimación de parámetros cuando 
los componentes del error se correla- 
cionan para un sistema de ecuaciones 
(Rose y Lynch, 2001); además permite 


RESULTADOS 

Modelos para volumen total y comercial 

En la tabla 2 se observan las esta- 
dísticas de ajuste y parámetros estimados 
para el modelo de volumen total sugerido 
por Schumacher y Hall (1933), con datos 
de remediciones en sitios permanentes 
de muestreo; el modelo presentó un 
ajuste significativo (p < 0.0001) y valores 
altos de R 2. 

Al ajustar simultáneamente los 
modelos de volumen comercial con el 
modelo de volumen total se obtuvieron 
los valores de los estimadores que se 
muestran en la tabla 3. El ajuste sugiere 
que el modelo de Fang y Bailey (1999) es 
el más adecuado para estimar el volumen 
comercial y total, con una mayor R2 y 
menores valores en la suma de 
cuadrados del error (sce), aunque el valor 
del cuadrado medio del error (cme) obte- 
nido en este modelo es igual al del 
modelo de Clutter (1980); sin embargo, 
se toma el modelo de Fang y Bailey por 
ser más sencillo en su estructura y con 
sólo cuatro parámetros a estimar. 


Tabla 2. Resultados del ajuste para el modelo de volumen total por mínimos cuadrados. 


Modelo R^ 

SCE 

CME 

Parámetros 

Valores 

estimados 

Error 

estándar 

Prob>T 





de tt, 





«0 

0,000032 

3,31 E -06 

<,0001 

1 0,9946 

0,00173 

0,000022 

a -¡ 

1,942372 

0,054 

<,0001 


«2 


0,911162 


0,058 


<,0001 
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Tabla 3. Resultados del ajuste para los modelos de volumen comercial por mínimos 

cuadrados. 


Modelo 


SCE 

CME 

Parámetros 

Valores 

estimados 

Error 

Estándar 

Prob>T 





ao 

0,00003 9 

2,49E-06 

<,0001 





ai 

2,20473 1 

0,0296 

<,0001 

2 

0,9814 

0,0244 

0,00002 5 

a2 

0,54628 8 

0,0302 

<,0001 





fio 

0,35699 

0,0223 

<,0001 





fh 

2,19903 7 

0,0304 

<,0001 





h 

-1,94498 

0,0331 

<,0001 





ao 

0,00003 6 

1,17E-06 

<,0001 

3 

0,9817 

0,0241 

0,00002 5 

ai 

1,94752 7 

0,0164 

<.0001 





a2 

0,82509 9 

0,0174 

<.0001 





fio 

0,00001 6 

1.13E-07 

<,0001 





ao 

0,00003 7 

1,09E-06 

<,0001 





ai 

1,97579 3 

0,0148 

<,0001 

4 

0,9809 

0,0251 

0,00002 6 

a2 

0,83256 5 

0,0157 

<,0001 





fio 

0,00001 9 

1,70E-07 

<,0001 





fh 

2,18589 1 

0,0303 

<,0001 


Durante el proceso de ajuste de los 
modelos para volumen total y comercial 
no se detectaron diferencias estadísticas 
significativas (p > 0.05) entre C. odorata y 
T. donnell-smithii, por lo que se tiene un 
solo modelo para ambas especies. 

Crecimiento en altura dominante 

El modelo Chapman-Richards 
presentó el mejor ajuste para predecir y 
explicar el crecimiento en altura de C. 
odorata y T. donnell-smithii {Tab\a 4), obte- 
niendo el mayor valor de Pseudo- R2, así 
como los valores más pequeños en la 
suma de cuadrados del error (sce) y en el 
cuadrado medio del error (cme). 

Al parámetro ói en el modelo (8) de 
Chapman-Richards le fue asignado un 


valor constante probando entre valores 
positivos menores de 1, y se observó 
que el valor de 0.08 fue el que mejor 
contribución tuvo al ajuste del modelo. 
Con el ajuste del modelo de Chapman- 
Richards se generaron familias de 
curvas anamórficas de is para 14, 18, 22 
y 26 m a una edad base de 10 años 
(Figura 2). 

El modelo de Chapman-Richards 
tiende a subestimar el crecimiento en 
sitios de menor calidad (is= 14) para 
edades tempranas y lo sobreestima en 
sitios de mejor calidad tanto para edades 
tempranas como posteriores (is= 26) 
(Figura 2). Sin embargo, estas curvas 
muestran que las especies estudiadas 
aún están en la fase acelerada de creci- 
miento en altura. 
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Tabla 4. Resumen del análisis de regresión efectuado con datos de altura-edad en 
Cedrela odorata y Tabebuia donnell-smithii para conocer el nivel de ajuste de los 

modelos empleados. 


Modelo 


SCE 

CME 

Parámetros 

Valores 

estimados de 6 ; 

Error 

Estándar 

Prob>T 


0,8707 

21,8289 

1,7463 

^0 

21,4092 

1,5001 

<,0001 

(6) 




Si 

3,923788 

0,3478 

<,0001 


0,919 

17,2047 

1,3764 

6 o 

29,6996 

2,3474 

<,0001 

(8) 




6 i 

0,08 

- 

- 





62 

1,241143 

0,09 

<,0001 



Edad en años 


Figura 2. Curvas anamórficas de índice de sitio para Cedrela odorata y Tabebuia 
donnell-smithii (edad base= 10 años) en plantaciones de enriquecimiento San José 

Chacalapa, Pochutia, Oaxaca. 
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La predicción del crecimiento en 
altura dominante no presentó diferencias 
estadísticas significativas (p > 0.05) entre 
las especies estudiadas, por lo que las 
curvas de índice de sitio pueden ser utili- 
zadas indistintamente para C, odorata y 
T. donnell-smithii. Sin embargo, gráfica- 
mente esta última especie muestra ligera- 
mente un mejor índice de sitio a edades 
tempranas (is= 26). 


Ajuste del SCRM 

Una vez definida la familia de curvas 
para altura dominante (A) se decidió 
realizar el ajuste simultáneo de ésta, así 
como del área basal (ab) y volumen {V) 
bajo la técnica sur (Tabla 5), la cual 
mostró un mejor ajuste en comparación 
con los mínimos cuadrados ordinarios; los 
modelos de proyección presentaron una 
alta R2, así como valores bajos de la sce 
y del CME respecto a los modelos de 
predicción (Tabla 5). 


Tabla 5. Parámetros estimados y bondad de ajuste bajo regresión aparentemente no 
correlacionada (sur), para los modelos de altura dominante (A), área basal (ab), y 

volumen total (V). 


Modelo 

Elemento 

G.L. 

model 0 

G.L. 

error 

SCE 

CME 


Re- 
ajustad a 

8 

Ai 

1,5 

12,5 

17,0802 

1,3664 

0,8961 

0,8920 

9 

A 2 

0,5 

13,5 

12,3438 

0,9144 

0,9173 

0,9204 

1 0 

ABi 

1 

1 3 

1 1 ,8424 

0,91 1 

0,8767 

0,8767 

1 1 

AB 2 

1 

1 3 

7,0121 

0,5394 

0,9484 

0,9484 

1 2 

Vi 

2 

1 2 

9,021 

0,7517 

0,9926 

0,9919 

1 3 

V 2 

1 

1 3 

7,0558 

0,5428 

0,9968 

0,9968 


Parámetro 

Valor 

Estimado 

SUR 

Error Estándar 
Aproximado 
SUR 

Valor de 
t SUR 

Valor de 
Rechazo SUR 

^0 

29,6468 7 

2,1303 

13,92 

<,0001 

Ó2 

1,22798 5 

0,0852 

14,41 

<,0001 

00 

5,23329 4 

0,5179 

10,11 

<,0001 

01 

0,98860 4 

0,0284 

34,76 

<,0001 

yo 

1 ,98647 1 

0,2116 

9,39 

<,0001 

yi 

1,105139 

0,3638 

3,04 

0,0103 

Z2 

1,29325 8 

0,0346 

37,38 

<,0001 






Tabla 6. Edad, altura dominante, área basal, volumen total, incremento medio anual (ima) e incremento corriente anual (ica) para cuatro 
índices de sitio (is) en plantaciones de Cedrela odorata y Tabebuia donnell-smithii en San José Chacalapa, Pochutia, Oaxaca. 
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La tabla 6 presenta la producción 
para las especies estudiadas a densi- 
dades iniciales variables entre los sitios 
de muestreo; en ella se representan, las 
alturas dominantes, así como el área 
basal y volumen por hectárea, además 
del incremento medio anual y corriente 
anual por hectárea por año a diferentes 
índices de sitio (Tabla 6). 

El crecimiento en volumen y área 
basal continúa de forma ascendente 
(Figura 3), situación similar muestran 
las curvas del incremento medio anual 
(ima) e incremento corriente anual (ica) 
en volumen (Figura 4). Es lógico pensar 
que debido a las edades jóvenes en las 
que se evaluaron las plantaciones, 
éstas todavía se encuentran en una 
etapa de crecimiento casi exponencial. 
Remediciones de estas parcelas sin 
duda arrojarán información muy valiosa 
sobre la duración de este comporta- 
miento. 

DISCUSIÓN Y CONCLUSIONES 

El modelo Schumacher y Hall 
(1933) empleado para definir el sistema 
de cubicación y obtener el volumen total 
por árbol en este estudio, fue utilizado 
por Reyes et al. (2006), para dos espe- 
cies de Eucalipto plantadas comercial- 
mente en Tehuantepec, Oaxaca, México, 
también por Ferreira et al. (2004) al 
emplear varios modelos matemáticos 
para obtener el volumen de Tabebuia 
cassinoides en Sao Paulo, Brasil y por 
Akindele y LeMay (2006), para un grupo 
de 77 especies de árboles tropicales, en 
Nigeria, aunque en especies y sitios dife- 
rentes, el modelo fue el que mejor se 
ajustó a los datos. 

Los resultados obtenidos en este 
trabajo difieren de los reportados por 
Bravo et al. (2007), quienes al desarro- 
llar un sistema de crecimiento y rendi- 


miento para caoba y cedro en planta- 
ciones de Tuxtepec, Oaxaca, emplearon 
los modelos de volumen total y comercial 
de Clutter (1980) y Piennar y Rheney 
(1988), así como el modelo de creci- 
miento en altura dominante de Schuma- 
cher. Mientras que en este trabajo estos 
modelos no presentaron los mejores 
indicadores de ajuste y estas diferencias 
posiblemente se deban al origen de los 
datos. 

Las plantaciones evaluadas en 
Tuxtepec, aunque de similar edad, no son 
de enriquecimiento de acahuales, sino 
plantaciones bajo silvicultura intensiva y 
bajo la mezcla cedro rojo-caoba. 

La representación gráfica del 
modelo Chapman-Richards en su versión 
anamórfica, para los datos en este 
estudio, difiere de lo reportado por García 
et al. (1998), quienes desarrollaron ecua- 
ciones de índice de sitio para Swietenia 
macrophylla King, en plantaciones en 
Quintana Roo, México, reportando que el 
modelo Chapman-Richards sobreestima 
el crecimiento para edades tempranas y 
lo subestima al final, esto es posible 
considerando las especies, sus hábitos 
de crecimiento y las edades de las planta- 
ciones, además esto es lo esperado dado 
que los árboles aún están en la fase 
acelerada de crecimiento. Sin embargo, 
en las mejores calidades de sitio se 
tiende a sobrestimar la altura; pero no es 
tan conservadora como el modelo de 
Schumacher, el cual tiende a subestimar 
la altura dominante en los mejores 
índices de sitio. 

El sistema de crecimiento y rendi- 
miento maderable desarrollado puede ser 
usado indistintamente para las dos espe- 
cies plantadas {Cedrela odorata, Tabe - 
bula donnell-smithii) y es una herramienta 
útil para la toma de decisiones en la 
planeación de su manejo futuro. 
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0 2 4 6 8 10 12 14 16 

Edad (años) 


Figura 3. Datos de área basal (A) y volumen (B), por índice de sitio, a una edad base 

de 10 años. 
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Edad (años) 


Figura 4. Curvas del incremento medio anual (ima) e incremento corriente anual (ica) 

por índice de sitio, a una edad base de 10 años. 


Por tratarse de una plantación joven, 
las especies estudiadas aun se encuen- 
tran en una fase acelerada de creci- 
miento, tanto en área basal como en 
volumen, situación que se confirma al 
observar la curva del incremento medio 
anual que todavía no llega a su punto de 
culminación. 

Para poder validar los estudios, a 
largo plazo, es conveniente dar segui- 
miento a los sitios de medición estable- 
cidos, por lo que los resultados presen- 
tados no deben extrapolarse más allá de 
los 13 años. Mediciones futuras permi- 
tirán establecer además patrones de 
mortalidad natural para el manejo 
adecuado de la densidad. 


El SCRM que se presenta es además 
auxiliar para los datos de inventarios 
provenientes de las plantaciones, ya que 
permite no sólo calcular el volumen por 
hectárea de manera exacta, sino proyec- 
tarlo con una confiabilidad aceptable. 
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